arXiv: 1505.008llvl [astro-ph.GA] 4 May 2015 


Annu. Rev. Astron. Astrophys. 2015 53 


Revisiting the Unified Model of 
Active Galactic Nuclei 

Hagai Netzer 

School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel, email:netzer@wise.tau.ac.il 


Key Words 

Active galaxies, unification, black holes, accretion disks, central torus, star formation, AGN 
surveys, interferometry, mergers 

Abstract 

This review describes recent developments related to the unified model of active galactic 
nuclei (AGN). It focuses on new ideas about the origin and properties of the central ob¬ 
scurer (torus), and the connection with its surrounding. The review does not address radio 
unification. AGN tori must be clumpy but the uncertainties about their properties are still 
large. Todays most promising models involve disk winds of various types and hydrody- 
namical simulations that link the large scale galactic disk to the inner accretion flow. IR 
studies greatly improved the understanding of the spectral energy distribution of AGNs 
but they are hindered by various selection effects. X-ray samples are more complete. A 
basic relationship which is still unexplained is the dependence of the torus covering factor 
on luminosity. There is also much confusion regarding “real type-II AGNs” that do not 
fit into a simple unification scheme. The most impressive recent results are due to IR in¬ 
terferometry, which is not in accord with most torus models, and the accurate mapping of 
central ionization cones. AGN unification may not apply to merging systems and is possibly 
restricted to secularly evolving galaxies. 


1 THE AGN FAMILY: CLASSIFICATION AND UNIFICATION 
1.1 AGN Classification 

An active galactic nucleus (AGN) is defined here as a galaxy containing a massive (> 10 s 
M@) accreting black hole (BH) with Eddington ratio exceeding the (somewhat arbitrary) 
limit of I/AGN/TBdd= 10~ 5 , where Lagn is the bolometric luminosity and L Edd = 1-5 x 
10 38 Mbh /Mq erg s ^ 1 is the Eddington luminosity for a solar composition gas. This defi¬ 
nition excludes the milky way Galaxy but includes a large number of low ionization nuclear 
emission-line regions (LINERs) and low ionization radio galaxies. Most AGNs include sev¬ 
eral of the following components: 

• A sub-pc rotational dominated accretion flow that is usually refers to as an accretion 
disk. Optically thick disks can be geometrically thin (“thin” accretion disk) or thick 
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(“slim” or “thick” accretion dwQ. Optically thin accretion disks or flows that are 
occasionally advection dominated. Such structures are referred to as radiati on ineffi¬ 
cient accretion flow (RIAF), or advection dominated accretion flow (ADAF) (INaravan 


2004 . 


• High density, dust-free gas clouds moving at roughly Keplerian velocities at a lumi¬ 
nosity dependent distance of 0.01-1 pc from the BH (the broad line region, BLR) 

• An axisymmetric dusty structure with luminosity dependent dimensions of 0.1-10 pc 
(the central “torus”). 

• Lower density, lower velocity ionized gas (narrow line region, NLR) extending from 
just outside the torus to hundreds and even thousand of parsecs along the general 
direction of the opening in the torus (“ionization cones”). Most of this gas contains 
dust except for very close in, in a region referred to as “the coronal line region”. 


• A very thin molecular maser disk with size similar to that of the torus. 

• A central radio jet occasionally associated with 7 -ray emission.. 


(Robson 19961. Peterson 1997 

Kroliki 1998. Blandford et al. 199d. Osterbrock & Ferland 

2006. Netzer 2013. Antonucci 

1993. Ho 2008. Heckman & Best 2014)1 and the reader is 


referred to these references for more information. 

“AGN unification” is the idea that the large diversity of observed AGN proper ties can be 


explained by a small number of physical p arameters. The old unification scheme ( Antonucci 
19931 . lUrrv fe Padovanilll995l. IUrrvll 200 3) is a courageous attempt to combine the ever in¬ 


creasing number of sub-groups into a general picture with two parameters: the torus incli¬ 
nation to the line of sight (LOS) and the source luminosity (“unification by inclination”). 
This scheme is perhaps the simplest possible way to characterize the known fact that the 
nuclear continuum and emission line radiation of AGNs can suffer wavelength dependent 
scattering, absorption and reflection on the way out. This can take place in the torus, in the 
disk of the host galaxy, in stellar and nuclear outflows, and inside the BLR itself. Recently 
it was suggested that AGNs can be sep arated in a different wa y into two major groups: 
“radiative mode” and “jet mode” AGNs ( Heckman fc Best 2014 1. Most of the energy out¬ 
put in radiative-mode AGNs is in the form of electromagnetic radiation, and is a direct 
result of matter accretion through a central optically thick accretion disk. This group is 
referred to in the literature as Seyfert galaxies or QSOs. About 10% of the sources in this 
group are radio-loud, showing a highly collimated, relativistic radio jet and, occasionally, a 
7 -ray jet. Radiative mode AGNs are very efficient accretors with I/AGN/I/Edd> 0.01. An 
alternative name, based on the level of ionization of the NLR gas, is “high ionization AGNs” 
(referred to also as “high excitation AGNs”). The prime energy output of jet-mode AGNs 
is bulk kinetic energy transported in two-sided jets. Their typical Eddington ratio is much 
smaller and the jets are most likely powered via a RIAF. The members of this group are 
low luminosity radio galaxies and LINERs, and an alternative spectroscopic definition is 
“low ionization” AGNs. In the local universe the mean BH mass of objects in the second 
group is larger. 

T his review concerns mos tly radiative mode AGNs and the reader is referred to Ho ( 2008h 
and lHeckman fc Best! (120141 1 for more information about radio mode sources. The exception 


1 The terminology here is not very clear and both names are used for systems that are not geometrically 
thin. In this article I use the term “slim disk”. 
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are some LINERs with total electromagnetic radiation that equals and even exceeds, the jet 
kinetic energy. I will not discuss radio-loud, radiative mode AGNs although most of their 
properties are indistinguishable from those of the radio-quiet sources. This excludes about 
10% of radiative mode AGNs. 

The main sub-groups addressed in this review are: 

Type-I AGNs: Sources showing broad (1000 — 20,000 km sec x ) permitted and semi- 
forbidden emission lines and a bright, non-stellar, central point source visible at all wave¬ 
lengths that are not contaminated by stellar light. Almost all low to intermediate luminosity 
type-I AGNs show strong, high ionization narrow emission lines many of which are forbidden 
lines. Most papers refer to these sources as Seyfert 1 galaxies, or QSOs. Narrow emission 
lines are missing from the spectrum of many high luminosity type-I AGNs. Many papers 
make a distinction between several sub-groups of type-I AGNs like type 1.5, type 1.8 or type 
1.9. This refers to the relative intensity of the broad and narrow components of the Balmer 
lines, type 1.9 being the ones with the stronger narrow component. Here I refer to all such 
objects as type-Ii where “i” stands for intermediate (in some papers “i” refers to broad 
Paschen lines). This has been a source of much confusion since part of the spectroscopic 
characteristics have little to do with the conditions near the BH (e.g. reddening in the host 
galaxy) and weak broad lines can also be the result of line and continuum variations. 
Type-II AGNs: Sources containing strong narrow (300 — 1000 km sec *) NIR-optical-UV 
emission lines that are somewhat broader than what is observed in emission line galaxies of 
similar types. The lines show clear indications of photoionization by a non-stellar source. 
The strongest lines are [Ohi] A5007, [Nn] A6584, [On] A3727 , [Oiv] A25.9/im, [Nev] A3426, 
CIV A1549 and the hydrogen Balmer and Lyman lines. Like type-I AGNs, they also show a 
point X-ray source. Advances in mid-IR (MIR) spectroscopy resulted in attempts to identify 
type-II AGNs by the equivalent width (EW) of the silicate 9.7/rm absorption feature in their 
spectrum. Such methods are not in very good agreement with optical classification. Type- 
II AGNs are further divided into two subgroups. The first contains hidden type-I sources 
with broad emission lines seen in polarized light. The second group is less accurately 
defined and is occasionally referred to as “real type-II” or “true type-IL’ AGNs. They 
show similar width and excitation narrow lines but no detectable broad lines and little X- 
ray absorption. Their mean luminosity is below the luminosity of the type-II objects with 
hidden broad lines. In the local Universe, such objects represent about 30% of all type-II 
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2002, 
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20031 . Shi et al. 2Q10l . Brightman fc Nandral 2011b , 


201 
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LINERs: AGNs that are characterised by their low ionization, narrow emission lines from 
gas which is ionized by a non-stellar source. Typical strong emission lines in this group 
are [N n] A6584, [N il] A6584 and [S il] A6731, and the Balmer lines. I will only address 
LINE Rs with EW(Hq)> 3A t o avoid the confusion with gas ionized by post-AGB stars 
fe.g. lCid Fernandes et al. 2011 1. LINERs can be divided into type-I (broad emission lines) 
and type-II (only narrow emission lines) LINERs. So me but not al l LINERs show point-like 
X-ray and UV sources and UV and X-ray varia tions jMaoz 20071. Hernandez-Garcia et al. 


2013 1. Detailed reviews of LINERs are given in IHoI ( 2008l l and lHeckman fe Bestl (l2014l l. 


Lineless AGNs: This subgroup consists of AGNs with extremely weak, sometimes com¬ 
pletely undetected emission lines. They are recognized by the presence of a non-stellar 
central source and, occasionally, continuum variabi lity. This gener al group includes two 
distinct sub-groups : low luminosity sources (e.g. Tnmrp_et_al, 200SJb and very high lumi¬ 
nosity sources (e.g. Shemmer et al. 20ld . Meusinger fe Balafkan 2014 1. 
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Classifying type-I AGNs is relatively easy except for cases where the BLR is heavily red¬ 
dened by foreground dust. Identifying and classifying narrow line AGNs is more challenging, 
especially when no hard X-ray data are available, partly because pure star forming (SF) 
galaxies (also referred to as HII galaxies), that are more numerous than AGNs, also show the 
same strong narrow emission li nes. An efficient way to separate the groups is to make use 


of “diagnostic diagrams” (see Baldwin 1 Phillips fc Terlevichl Il98ll : lYeilleux fc Osterbrock 


1987 : Kauffmann et al. 20031 : KewieT^Takl 2006 ) that compare various line ratios that dif¬ 


fer according to the nature of the ionizing continuum (stellar or non-stellar). The most 
useful line ratios are [O hi] A5007/H/3, [N n] A6584/Ha, and [O i] A6300 /Ha. Diagnostic di¬ 
agrams show several well defined regions, one including high ionization AGNs, one including 
mostly LINERs, and one occupied by SF galaxies. The region between the pure SF systems 
and AGNs is occupied by “composite sources”, with intermediate spectral properties. 

Most AGNs known to date were found in large surveys, like the Sloan Digital Sky Survey 
(SDSS), through spectroscopic observations with relatively larg e apertures (e.g. 3 arc-se c 
in the SDSS). The host galaxies of most AGNs are SF galaxies ( Rosario et al. 2012, 20131 '). 
We can estimate the total emitted SF luminosity, Lsf, and distinguish “AGN dominated” 
(Lagn>Lsf) from “SF dominated” (Lsf>Lagn) systems. For example, in the spectroscopic 
3 arc-sec SDSS sample, such an comparison shows that many composite sources are SF 
dominated galaxies. Difficulties in identifying AGNs in composite sources can affect the 
census of local type-II AGNs, their luminosity and mass functions (LFs, MFs). The census 
of type-I AGNs is based on other methods and is not subjected to such uncertainties beyond 
a redshift of about 0.1. 


Redshift 

SF galaxies 

LINERs 

Type-II AGNs 

Composite sources 

0.05-0.1 

0.755 

0.035 

0.04 

0.17 

0.15-0.2 

0.605 

0.035 

0.11 

0.25 


Table 1: The fraction of strong emission line sources, at two redshift intervals, as derived 
from diagnostic diagrams of SDSS objects from the data release 7 of the catalog. The total 
number of sources in the 0.05-0.1 redshift interval is 108,353 and the number in the 0.15-0.2 
interval, 17,572. The selection includes all emission line objects where the lines that were 
used for the classification have signal-to-noise ratio larger than 3. This introduces a large 
bias against sources at higher redshifts, especially weak line LINERs (courtesy of David 
Rosario). 


1.2 The Observed Signature of a Central Obscurer 

1.2.1 O ptical -UV Absorption and Reflection The original unification scheme 
( Antonucci 19931 . and references therein) included a central axisymmetric dusty structure 
with a column density which is large enough to completely obscure the central source in 
some directions. This is usually referred to as the central “torus”; a somewhat ambiguous 
name which will, nevertheless, be retained here. The simplest torus is made of a smooth 
matter distribution. More elaborated structures, made of clumps and inter-clump material, 
are preferred by today’s observations. The gas at the inner radius of the torus is ionized by 
the central source. Deeper in, the torus contains dusty molecular gas. 

It is interesting to study the observational consequences of placing such a simple geometri¬ 
cal structure around the central BH and accretion disk. The first prediction is a dependence 
of the obscuring column on inclination to the LOS, where inclination is measured relative to 


4 Netzer 







































the z axis of the system (the direction perpendicular to the central disk). Small inclinations, 
close to face-on, allow a clear view of the central source and high inclinations result in large 
obscuration at UV-optical wavelengths. 

The second prediction is based on the fact that electromagnetic radiation reflected off 
dust grains and free electrons is polarized, with angle and percentage polarization wh ich 
depend on geometry and wavelength. The seminal work of lAntonucci fe Milled (119851 1 is 
among the first of many papers showing how to detect hidden broad emission lines using 
spectropolarimetry and how to infer the BLR ki nematics a nd overall geometry from such 
observations. Later and deeper polarimetry (e.g. lTra n l2003T ) show a large fraction of type- 
IIs with no hint for polarized broad lines. This led to the suggestion that objects lacking 
polarized broad lines are unobscured but void of BLR gas, like the real type-II AGNs fS ll.lt 
that show no detectable broad lines even when not obscured.. 


1.2.2 X-ray Absorption Reflection and Emission The next unavoidable conse¬ 
quence is a column density dependent X-ray absorption. The assumption used here is of an 
isotropically emitting central X-ray source, possibly a hot corona over the central part of 
the accretion disk. For a Compton thin torus, with Nh < 1-25 x 10 24 cm" 2 , the obscuring 
column density can be measured by modeling the observed spectrum provided the intrinsic 
spectral energy distribution (SED) and the gas metallicity are known (for solar metallicity 
gas, most of the X-ray absorption beyond 0.3 keV is due to metals). 

There are several X-ray spectral features that can help diagnose the torus properties. 
The first is the “Compton hump”, an excess radiation at around 20 keV due to the down 
scattering of hard X-ray photons. Almost all type-I AGNs show Compton humps in their 
spectra (the terms “scattering” and “reflection” are used quite l oosely i n the literature, all 
referring to Compton scattering followed by absorption, see e.g. Netzer (1993!)). 

The second feature is due to reflection off electrons in highly ionized gas. Here the 
reflected radiation is a mirror image of the intrinsic continuum with intensity that, for all 
practical purposes, depends only on the geometrical covering factor (/c=fl/47r where f l is 
the solid angle subtended by the torus). This fraction hardly ever exceeds a few percent. 

The third spectral feature is the K a iron line. This emission feature is a collection 
of fluorescence (Auger) K a transitions and recombination of H-like and He-like iron at 
energies between 6.4 and 7 keV, depending on the specific ion. The lines are broad and of 
high ionization if originating from the surface of the central accretion disk, or narrow and 
occasionally of low ionization if due to X-ray illumination of the cold gas in the torus. In 
case where the X-ray continuum is unobscured, the EW of the narrow Ka line depends on 
the number of absorbed E> 7.1 keV photons. For a typical X-ray SED and solar metallicity 
gas, this is estimated to be: EW(Ka)«/cxATr/10 21 eV, where complications like multiple 
scattering have been ignored (note that Nh is not necessarily the LOS absorbing column, see 
§3.3) For a torus with a covering factor of 0.5 and column density of 10 23 cm " 2 , EW(Kq' 


50 — 1 00 eV. Indeed, m any type-I AGNs show narrow Ka lines with this EW (|Bianchi et al 


20081 . iRicci et alJboij . and references therein). An even clearer signature of a central torus 
is the very large EW(Ka) in Compton thick sources. In this case, much of the Ka line 
reaches the observer without being absorbed but most of the central continuum radiation 
is blocked by the torus. The observer can only see the much weaker reflected continuum 
which raises the observed EW dramatically and it can reach 1-2 keV. 


1.2.3 Ionization Cones The unobscured ionizing radiation of the central source 
reaches the gas in the galaxy causing heating and ionization. This results in the emis- 
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sion of various narrow lines. A torus like geometry will result in two “ionization cones” (or 
bicone) Depending on the geometry and inclination to the line of sight, one of the cones 
may not be observed. 


1.2.4 The Dusty Boundary of the Broad Line Region The dust sublimation 
temperature, T su b, is the maximum temperature attained by a dust grain before it evap¬ 
orates. This depends on grain size and composition and the local flux. For ISM type 
graphite grains, the maximum sublimation temperature is T su b — 1800 K and for silicate 
grains, T su b — 1400 K. These numbers can be used to derive a mean sublimation radius, 
(Rsub), which is the minimum radius where a grain of a certain size can survive the lo¬ 
cal radiation field. Averaging over ISM-type grain sizes, one obt ains the following mean 
sublimation distances for graphite (C) and silicate (Si) grains fe.g. lBarvainisHl987l) : 


(Rsub.c) ^ 0.5L ^ 2 


1800 


T s , 


f(0) pc; (Rsub,si) — l- 3 I/ 4 g 2 


1500 




/W pc, (1) 


where Z,46 =Lagn/10 46 erg s~ and f(6) is an angular dependent ter m that allows for 
anisotropy of the central source radiation ( Netzer fe Trakhtenbrot 20141) . For an isotropic 
source f (9) = 1 and for a thin disk with electron scattering atmosphere, f(9) = [cos 6(1 + 2 cos 9) /3] 
(Netzei 20131 ). The above numbers are computed for a mean grain size of a — 0.05/xm and 
the distance depends roughly on a -2 . 

Dust can compete effectively with the ionization of the gas because of its large absorp¬ 
tion cross sections at all wavelengths longer than about 0.02 /jrri. For highly ionized gas, 

A r d„st /A’ h o oc IVgas/Ajjo oc Uu, where U h is the ionization parameter defined over the hy¬ 
drogen Lyman continuum ( Netzei 2013l h At high levels of ionization, the dust absorbs a 
larger fraction of the photons capable of ionizing hydrogen and helium. This reduces the 
size of the ionized part of the cloud and lowers the intensity of most emission lines. Emission 
line photons produced far from the surface of an ionized, dusty gas cloud, suffer additional 
attenuation on their way out due to absorption by the dust. This is most noticeable in lines 
with high optical depth where the large number of scatterings of the line photons increase 
considerably the path length. Thus, a clear prediction of the torus scenario is that the high 
velocity gas near the BH will be divided into two dist inct regions, dust-fre e gas with strong 
line emission and dusty gas with weak line emission ( Netzer fc Laor 1993). 


1.2.5 Infrared Emission A clear prediction of any torus model is a broad IR SED 
with a total luminosity that depends almost entirely on the fraction of the central source 
radiation absorbed by dust in the torus (a combination of optical depth and covering factor). 
The radiation may be emitted isotropically, or anisotropically, depending on the optical 
depth and exact geometry. The SED is predicted to show several silicate features centered 
around 9.7 and 18/tm, either in absorption or emission, depending on the exact geometry 
and LOS optical depth. The emitted radiation of basically all AGNs is varying in time. This 
will result in a time-dependent incident flux on the inner torus walls, and time dependent 
NIR and MIR emission by the dust. 


1.3 Unification Schemes 

The term “unification scheme” refers usually to two different categories: IR-optical-UV- 
X-ray unification, and radio unification. The first involves a central dusty torus and can 
explain the major differences between type-I and type-II AGNs with a surprisingly small 
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number of assumptions. This was first suggested in the late 1980s. The second combines 
the torus with a relativistic jet observed in some 10% of all high ionization AGNs. The jet 
is launched in the vicinity of the BH and is aligned with the symmetry axis of the system. 
This can be used to make specific predictions about the observed properties of compact and 
extended ra dio-loud AG Ns. Radio unification is not covered in this review and the reader 
is referred to Urrv 1 20031 ) and Tadhunter ( 20081 ) for references on this topic. It is also worth 
noting that LINERs were not included in the first unification schemes. 

Recent observations of large AGN samples, and more detailed information about specific 
sources, lead to the conclusion that the old unification scheme requires three major mod¬ 
ifications. The first is related to the physics and structure of the central obscurer. The 
second is the realization that there are real differences between various sub-groups, beyond 
the dependence on luminosity and torus inclination. This can be due to the nature of the 
central power-house (e.g. a RIAF instead of high efficiency accretion flow), and the fact 
that some AGNs lack one or more of the main components, e.g. a BLR or a NLR. The third 
is related to BH and galaxy evolution. The present review focuses on these modifications. 


1.4 The Era of Selection Effects 


Understanding the cosmological evolution of super-massive BHs is key to the understanding 
of galaxy evolution This requires reliable luminosity and mass functions which depend on 
AGN classification. While large surveys like SDSS contain hundreds of thousands of AGNs 
up to z ~ 7, this survey, and basically all others, are hindered by various selection effects. 
Since the topic of this review is AGN unification, I will only mention here selection effects 
that can bias AGN classification. 

The use of diagnostic diagrams in spectroscopic surveys, like SDSS, to obtain complete 
AGN samples is limited to low redshift. Beyond z ~ 0.2, type-II AGNs with weaker emission 
lines drop from the sample because of poor signal to noise (S/N) and contamination by 
stellar absorption lines (weak line LINERs are already missing from such samples at redshift 
as small as 0.1). At z > 0.4, the Ha line, which is crucial for the classification, is shifted 
beyond the o bservable range. A n alternative method to detect luminous type-II AGNs up 
to z ~ 0.83 I Reves et al. 20081 . and refferences therein) is to search for strong [O ill] lines 
relative to H/3 which is required not to show a broad component. Unfortunately, large IR 
spectroscopic surveys are not yet available, and alternative diagnostics, based on emission 
lines at shorter wavelengths, are not as reliable. Broad line AGNs are detected by their 
continuum shape and their broad emission lines. At z < 0.1, the least luminous type-Is 
drop from such samples due to strong stellar continuum and large amplitude continuum 
variations. 

Reddening and obscuration is another reason for losing sources from large samples. At 
small redshift, this affects type-IIs more than type-Is. IR-selected samples are less prone 
to reddening but are biased towards sources with larger covering factor and against large 
inclination angles (see (j2j. They are also contaminated by emissio n from SF region s in 


the galaxy. Mid-IR photometric surveys, like those based on Spitzer (IStern et al. 


very efficient in differentiating AGNs from non-active galaxies. However, as will be shown 
later, the IR-SED is not uniform enough to unambiguously differentiate between AGN 
types and type-II classification is most problematic. Host galaxy contamination affects the 
measuring of low luminosity AGNs over the entire rest-frame range of 0.3-2 pm. Selecting 
and classifying AGNs by their hard X-ray continuum is, arguably, the most reliable way at 
almost all redshifts. The completeness of this method depends on the energy and flux limit 
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of the specific sample and various sophisticated methods have been developed to allow for 
such biases. A major difficulty is to detect highly obscured, large column density X-ray 
sources at high redshift. 

A co mprehensive comparison of the efficiency of various methods to discover AGNs is 
given in 


Merloni fe Heinz ( 20131 ). According to this work (see their Figure 16 and related 


to 90% at Lagn=10 42 erg sec 1 . 
to 70% at Lagn=10 42 erg sec -1 . 


explanation), the missing fractions of AGNs in various samples depend on the bolometric 
luminosity and are roughly: Optical (B-band) samples, from 70% at Lagn=10 47 erg sec -1 

MIR (15/rm) samples, from 50% at Lagn=10 47 erg sec -1 
. Hard X-ray (2-10 keV) samples, from 0% at Lagn=10 47 
erg sec -1 to 40% at Lagn=10 42 erg sec -1 . 

The following sections review todays torus models (§2), and torus observations (§3). 
This is followed by a description of additional, smaller scale and larger scale obscurers (§4), 
ionization cones (§5) and other galactic components connected to the torus (§6). Finally, in 
§7, I explain why AGN unification may not apply to major mergers. §8 is a short summary 
of the main points of this review. 


2 DISK AND TORUS MODELS 
2.1 Central Disks and Central Tori 

How large and how massive is the central structure? How is it connected to the BLR and 
accretion disk on its inner side, and to the ISM of the host galaxy on the outer side? Is it 
a torus, a disk, or some other shape? Is such a structure typical of all large bulge galaxies 
or is BH accretion essential to form and maintain it? And most important, how is this 
structure formed in active galaxies? The global theoretical framework is presented in this 
section and the observational consequences in §3. 

The terminology used in this area of research requires some clarification since the terms 
“central disk”, “central torus” and “central accretion disk” have been used quite loosely. In 
this review, the term inner accretion disk refers to the accretion flow inside the sublimation 
radius, where the accreted gas is losing angular momentum through molecular, magnetic or 
other sources of viscosity. Many AGNs must contain a much larger disk-like structure, 1-100 
pc in size, that connects the inner disk to the host galaxy. The inner part of this structure 
can have a thick geometry in which case it becomes the inner torus. The entire structure, 
that can extend beyond the torus, is referred to here as the “central disk” or the U Q disk” 
(after the Toomre stability parameter Q, see below). It contains a large amount of molecular 
gas that, given a large enough column density, undergoes SF associated with super-novae 
(SN) explosions, the growth of a nuclear star cluster, etc. Some of these processes result 
in local effective viscosity that can make this entire structure a large accretion disk, where 
gas loses its angular momentum and drifts into the center. The angular momentum of the 
large disk is not necessarily aligned with the angular momentum of the inner disk. There 
is also no known reason why such a large disk will maintain its properties throughout the 
accretion episode. 


2.2 General Models 


Two general types of central disks and central tori have been considered. The first are theo¬ 
retic al models (e.g. Krolik_&Begehnan|l98S JCollhi&_Zahn|l999t lEmmeri ng. Blandford fe Sh l osmanl 

J: 


1992: Schartmann et al . 2005:lcollin fc Zahnl2008ljBeckert fc Duschll2004l;lHonig fe Beckert 
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2007; Krolild2007 ; Ehtzurj^_SI^^ Vollmer. Beckert fc Daviesl2008l : Wada. Panadopoulos fc Spaan; 

ta r 


2009: IWada 12012 : ISchartmann et al.ll201C ) that were developed to explore specific aspects 
like the formation, stability, shape and evolution of the central structure. The models 
attempt to answer two fundamental questions: what is the mechanism that removes an¬ 
gular momentum from galactic infalling gas thus allowing the feeding of the central BH, 
and what are the ways to maintain a thick gaseous structure for a long period of time. 
Important processes considered in this context include blob motions along magnetic field 
lines, colliding magnetic clumps, radiation pressure supported dusty clumps, SF driven 
turbulence, SN explosions, AGN and stellar feedback via radiation pressure and winds, 
and more. Several analytic and semi-analytic models have been proposed but the hydro- 
dynamical simulations seem to provide more realistic ways to model this complex situa¬ 
tion. The second type is phenomenological, ad-hoc models that make assumptions about 
the size, composition and geometry of the torus in an attempt to explain the observa- 
tions, mostly the spatially resolved images and the NIR-MIR SEP ( e.g. Pier_&^JCrplik 


1992 


Fritz 1 jYanceschini fe Hatziminaoglou 20061 : Nenkova et al. 2008a : Honig et al.ll2010l : 
Stalev ski et alj 20121 '!. Most such models neglect the connection between the torus and its 


surrounding. 

A natural way to describe the accretion problem is to start from the inner disk, assumed 
in most models to be geometrically thin and optically thick, where the local viscosity allows 
radial drift of material all the way to the innermost stable circular orbit (ISCO). A natural 
outer boundary of the accretion disk is its self gravity radius, Rsg'■ the location where local 
gravity exceeds the vertical component of the central BH gravity and the disk becomes 
unstable (Q = 1). This radius is given by: 


Rsg ~ 1680Mg 2/9 a 2/9 [LAGN/L E ddf /9 [n/0.1}- 4/9 R g 


( 2 ) 


where M 9 =Mbh/ 10 9 Mq, a is the viscosity parameter which is of order 0.01-0.1, 77 is the 
mass -to-radiation conversion efficiency, and R g is, again, the gravitational radius (ILaor fe Netzei 
1989lf . For a stationary Mg = 1 BH and a thin accretion disks radiating at Z/AGN/TEdd~ 0.1, 
Rsg — 0.04 pc. 

In the absence of additional local heating sources, Q drops below unity beyond the self¬ 
gravity radius, and the disk is fragmented into clouds that move in the same general plane. 

The gap that is formed prevents t he inflow of additional, more distant material to the 
center. According to some models ( Collin fe Zahn 19991 . 20081 . Duschl fe Britsch 20061 ') . a 
marginally stable situation with Q ~ 1 can be maintained over a large distance beyond Rsg 
due to collisions between clouds, and SF and SN explosions that produce turbulence and 
local viscosity. This allows accretion from much larger radii. For very massive BHs, such a 
marginally unstable disk can extend to ~ 100 pc and beyond. A challenge of models of this 
type is to combine all these processes into a coherent SF history which allows to follow the 
evolution of the nuclear star cluster. 

Another comprehensive semi-analytic model by ( Vollmer. Beckert fc Davies! 2008 1. con¬ 
siders gas infall from the galaxy which results in the formation of a clumpy central disk. 
This leads to the formation of a torus through three phases of evolution: 1) An initial 
short and massive gas infall leading to the formation of a turbulent and stellar wind-driven 
Q ~ 1 disk. This is followed by energetic SF processes. 2) Once supernovae explode, the 
inter-cloud medium is removed, leaving a massive, geometrically thick, collisional disk with 
a decreasing, but still high-mass accretion rate. 3) When the mass accretion rate has sig- 
nicantly decreased, the torus becomes t hin and transpa rent, similar to the circumnuclear 
disk in the center of the Milky Way (e.g. Alia et al. 2013l l. Numerical hydrodynamical sim- 
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ulations by Schartmann et al. ( 2010 1 follow the final stages of the model focusing on BH 
feeding during phases 2 and 3. The suggestion is that low velocity stellar winds, during the 
asymptotic giant branch (AGB) phase, starting some 50 Myr after the onset of SF, contain 
enough gas to feed a small central BH for a period of about 10 8 yr. The origin of such 
winds is a quasi-spherical star cluster with little angular momentum. The winds collide, 
form a central thin disk, and the gas is accreted onto the BH via an unspecified viscosity 
mechanism. Such models are currently limited to small, few x 10 7 M@ BHs. 

A gaseous torus cannot retain a large height/radius ratio with pure rotation. Star clus¬ 
ters can retain such a shape but not the gas between the stars since cloud collision will 
drive them to the plane on a few dynamical time scales. Maintaining a thick gaseous struc¬ 
ture requires high velocity turbulent motion, or outflow motion, and see ms to wor k better 
in a clumpy medium where there are frequent cloud-cloud collisions ( Beckert fc Dusch] 
2004). Other proposed m echanisms to maintain a thick structure involve UV, optical and 


IR rad iation pressure (e.g. Honig fc Beckert] 2007. Krolikl 2007. Czerny fc Hr^niewicz 


Wada 


2011 


2012|), magnetic winds (Emmering. Blandford fc Shlosman 19921 : Konigl fc Kartie 


1994 : Elitzur fc Shlosman 2006l l. and SF in the inflowing gas. UV-optical radiation by 
the central accretion disk, and locally (inside the torus) emitted IR radiation, can be com¬ 
bined into a plausible scheme that helps to solve the scale height problem. Here clumps are 
being pushed in the radial directio n by t he central so urce radiation in a way that depends 
on I/AGN/TEdd ( Honig fc Beckert 2007 . Wada 20121 4 . The central source heats the dust 
in the clumps which in turn emits IR radiation. The local IR flux diffuses outward being 
absorbed and re-emitted by neighboring clumps. This results in a radiation pressure forc e 
that balances or overcomes the z-component of the central source gravity (IKrolikl 2007). 

In this picture, an active AGN is necessary to both produce the toroidal structure and 

maintain its thick structure. _ _ 

The disk wind scenario propo sed in Emmering, Blandford fc Shlosman ]l992|) a nd dis- _ 

cussed in numerous other paper s ( Konigl fc Kaitjel l994l : Kartie fc Koniglll99f)l : Bottorff. Korista fc Shlosman 


2000l: Elitzur fc Shlosman! 200fil: EfitzurfcHof2009[ ) consists of a continuous large disk which 

is fed by cold material from the host galaxy. Close to the BH, it develops a magnetocentrifu¬ 
gal wind that supports rising clumps. In this scenario, the torus is merely a region in the 
wind that happens to provide a toroidal-clumpy structure. Clumps inside the sublimation 
radius are dust-free and can be viewed as BLR clouds. At larger distances, the clumps are 
dusty and can be viewed as part of the torus. The exact column density of the clumps differ 
from one model to the next. It may depend on the magnetic field properties, the ratio of 
the external accretion rate to the mass outflow rate from the disk, and more. The clumps 
that are most relevant to the observations are those capable of emitting strong lines, with 
a column density of at least 10 22 cm~ 2 . As the clumps rise from the disk, their column 
density drops dramatically and they disappear from view. The mass outflow rate in several 
wind models is somewhat smaller than the mass accretion rate through the disk. 

A recent development is the understanding that dust can survive in the atmosphere of 
the inner accretion disk at all distances where the disk eff ective temperature, T e //, is below 
the sublimation temperature ( Czerny fc Hrvniewic dEoil ). Simple estimates show that at a 
distance of Rblr (the emissivity weighted radius of the BLR), T e // ~ 1000K, roughly the 
dust sublimation temperature. This raises the possibility that, starting from this distance, 
the local radiation pressure acting on the local grains, is effective enough to raise dusty 
clumps from the disk surface. As the clumps rise, they are no longer shielded from the 
central radiation field which is more intense than the local field, they are being pushed in 
a radial direction, and the dust evaporates. This reduces the push in the vertical direction 
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Figure 1: 

Relative line emissivities per unit covering factor for twice solar metallicity BLR gas as a function 
of distance from the central BH. For the AGN luminosity considered here (7 x 10 46 erg s —1 ), the 
graphite sublimation radius is at 1.32pc and the silicate sublimation radius at 3.5pc. These 
distances are clearly visible due to the large drop in the intensity of all lines resulting from the 
absorption of the ionizing radiation by the grains, and the selective depletion of most metals. 

No te that some l in e emi ssion is still present well within the dusty torus (for details of the model 
see iMor &; Netzerl li2012T )l. 


and the clumps fall back onto the surface of the disk. Further away, the grains in the lifted 
material can survive the central radiation field and they form a clumpy thick structure. This 
“failed wind” scenario is claimed to explain several BLR observations. Its main uncertainties 
are related to the nature of accretion disks (thin, thick, etc), the variability of the central 
source radiation, and most importantly, the question of whether the local disk radiation 
pressure can support very large column density dusty clumps. 

Disk outflow models, magnetocentrifugal wind or radiation pressure driven wind, result 
in one continuous structure whose geometry is dictated by the global accretion process, and 
whose division into dusty and dust-free regions depends entirely on the central radiation 
field. A major distinction between the dust-free clouds closer in (BLR clouds) and the 
dusty ones further away, is the efficiency of line formation. This can be calculated by 
photoionization models that provide the predicted emissivity of a specific broad emission 
line as a function of distance from the central ionizing source. One example is shown in figure 
[T| The parameters of the dust-free clouds in this model were chosen to reproduces a typical 
broad line spectrum with the correct Rblr- The diagram shows the dramatic decrease 
in line emissivity at the graphite and silicate sublimation radii due to the appearance of 
grains which absorb most of the ionizing flux. The details of the line emissivity depend on 
the run of ionization parameter but the global behaviour must be typical. Note that line 
emissivities do not drop to zero and there is additional emission from the dusty part, inside 
the torus. 

Perhaps the most realistic models of today are numerical hydrodynamical simulations 
that include most but not all of the processes discussed above. T he t hree-d i mensi onal 
hydrodynamical simulations of Wada. Papadopoulos fe Spaansl ( 2009l 'l and Wada ( 2012 ) are 
perhaps the most detailed of their kind. They treat the ISM around a BH with ~ 10 7 M@ 
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Figure 2: 

Left: Face-on and edge-on density maps of a central toroidal structure resulting from gas infall 
towards a 1 (V MBH in a center of a spiral galaxy as calculated by 

IWada. Papadopoulos fe Spaanj (12009! ). Physical processes that are included in the simulations are 
SF, SN explosions and stellar feedb ack (Courtesy o f Keiic hi Wada). Right: Edge-on view of the 
2-phase phenomenological model of lStalevski et all d2012f ) showing high condensation clumps in 
yellow-orange colors and lower density inter-clump dust in blue (Courtesy of Marko Stalevski) 


and assume SF and viscosity which increases during the SN phase. The 2009 work considers 
only stellar feedback. It tracks atomic and molecular hydrogen with a spatial resolution of 
0.125 pc. In a quasi-equilibrium state, the gas forms an inhomogeneous disk with a typical 
diameter of about 30 pc. The velocity field of the disk is highly turbulent in the torus 
region and the average SN rate is small but enough to energize the thick structure. The 
computed gas column densities are consiste nt with tho se derived by X-ray observations. 
A more advanced version of the calculation ( Wada 20121 ') includes AGN radiation pressure 
(AGN feedback). This drives a “fountain” and the interaction between the non-steady 
outflow and inflow results in the formation of a geometrically thick turbulent torus. Face 
on and edge on views of the density in one such model are shown in figure [5] 

Present day numerical calculations have limited scope, in particular they treat low mass 
BHs and low accretion rate systems. It is not at all clear that such scenarios are applicable 
to AGNs with BH mass of order 10 9 Mq and accretion rates of order lMg/yr. At high 
accretion rates, the local stellar population may be different, the SF events more violent, 
and the stellar population in the central region different from the one observed in nearby 
systems.. 

It is important to note two other consequences of all torus models considered here. The 
first is the connection of the torus typical dimensions with the gravitational potential of the 
galaxy which can be quantified by considering the “sphere of influence” of the central BH. 
This sphere determines the relative importance of the BH gravity and the bulge gravity, on 
material in the center. It is defined as 


I?BH, S ph = - 10.7 


Mb 


r _^_r 

L 200km s 1 J 


10 8 Mq L 200km s 


pc, 


(3) 


where < 7 * is the stellar velocity dispersion in the bulge. For I/AGN/TEdd=0.1, this corre¬ 
sponds to few x 10 6 (<r*/200kms _1 ) _2 i? g . For nearby AGNs, I?BH,sph ~ 10 pc, similar 
and even smaller than the predicted outer dimension of the torus. Thus, the gravitational 
influence of the stars in the bulge is comparable to that of the central BH and must be 
included in the calculations. 

The second consequence is related to broad line observations. All the models considered 
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here require continuous flow of gas all the way from several pc to the ISCO. This would 
mean blocking the view of the part of the BLR behind the disk in type-I AGNs that are 
observed close to face-on. The observational implications of this scenario have not been 
studied so far. The masses in the torus and the central accretion disk are both of order 1% 
of the BH mass, which, given an accretion rate that produce I/AGN/LEdd=0.1, could supply 
enough gas for few x 10 6 yr. This raises the possibility that a long accretion episode of 
~ 10 s yr is broken into several shorter episodes during which the inner disk is disconnected 
from the outer flow. 


An inter esting variant of th e general torus formation scenario, proposed by IWang et al 


( 20ld l and Wang et al. ( 2011bl i , moves the SF activity to inside the dust sublimation radius. 
This model is not concerned with the issues of torus formation and stability. It merely 
assumes the existence of a central torus which serves as a source of mass supply to the 
innermost region. The model includes a thin, central accretion disk around the central BH, 
and fragmented clouds between Rsg and the inner walls of the torus. SF activity and SN 
explosions take place inside the torus walls and the duration of one accretion episode is 
estimated to be few x 10 6 yr. This would result in ~ 10 6 stars that may eventually settle 
into a central star cluster. 


2.3 Phenomenological torus models 


2.3.1 Continuous and Clumpy Tori Phenomenological torus models is the name 
chosen here for ad hoc models that are constructed to enable a more accurate treat¬ 
ment of the radiative transfer, and better modeling of the emitted SED and the images 
of the torus at different wavelengths. There are three generic type s of phenomenologi¬ 
cal torus models: continuous (or smooth) gas distribution models, (e.g. iPier fc Kroliklll992 


cal torus models: continuous (or smooth) gas distribution models, (e.g. ricr &z Kroiik 1992; 

Crmato J _^aues£_fc_jYajir^schinj|l997j^ Efatathiou fc Rowan-Robinsonll995 : Dullemon^&jai^^emnml 


200d), clumpy torus models, (e.g. DuflemoncTfc'Tair^emme] 


2005 


2005 


Nenkova et al 


12008a. b. Honig 


and continuous) models (IStalevski et al 


eta 


J201C ), and composite (a comb ination of clumpy 
20121 ). The Feltre et al.| (2012 ) wor k is an illu¬ 


minating comparison of the first two types and the work of Lira et alT ( 2013h provides 
comparison of two clumpy torus models based on fitting the spectra of type-II AGNs. 

All phenomenological models are axisymmetric and the torus inner walls are assumed to 
be at the dust sublimation radius. Unfortunately, the definition of this radius can differ 
between models because of the different types of grains, and the exact way used to estimate 


the bo lometric luminosity of the central sources. Some of the models, e.g. IStalevski et al. 
(2012), take into account the anisotropy of the central source radiation. This adds an 


angular term to the torus equations (/(#) in eqn. [lj to describe the shrinking sublimation 
radius at large polar angles. In this case, the torus equatorial plane reaches all the way to 
the outer boundary of the central accretion disk. The effect is more noticeable in the case of 
slim disks, where the angular dependence of the central source radiati on is more noticeable 
( Kawaguchi fe Mori 2010l . Kawaguchi fe Mori 2011 . Wang et al. 2014 ). 

The important parameters of the continuous phenomenological models are the inner and 
outer radii of the torus, the density profile and the opening angle. Additional parameters of 
clumpy tori models are the column density of individual clumps, the radial distribution of 
the clumps (this replaces the run of gas density in the smooth model), the filling factor and 
density distribution of the clumps. Present day phenomenological models do not yet include 
different types of grains in different parts of the torus although t his has been considered i n 
simpler (from the point of view of radiative transfer) models like Schartmann et al. ( 20081 ). 
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Radiative transfer is relatively easy to handle in continuous dusty tori but this is not the 
case in clumpy tori. Different st udies adopt dif f erent a pproaches to solving this problem. 
The analytical ID calculations of lNenkova et ali (|2008a ) are limited in their ability to treat 
complicated geometries but can cover a large parameter space wit h different c l ump sizes, 
torus opening angle a nd clump distribution. The 3D calculations of Honig et al. ( 20101 ) and 
Stalevski et al. ( 20121 '! are example of models using a Monte Carlo (MC) transfer technique, 


combined with ray-tracing, that allows more accurate solutions and, in principle, can deal 
with different grain size and composition. However, the calculations are computationally 
intensive and the parameter space explored more limited. 

Continuous dust distribution models are characterized by a monotonically decreasing 
dust temperature as a function of distance from the BH due to the F oc r~ 2 term and the 
attenuation of the incident radiation by the dust. In such models, every radius corresponds 
to a certain temperature. One prediction of continuous models is strong silicate 9.7/rm 
absorption in cases where the total optical depth at this wavelength exceeds unity. This is 
simply the result of the absorption of the flux emitted at this wavelength by the inner torus 
walls by further away material. This is pred icted to be most notic eable in type-II sources 
that are seen at large inclination angles (e.g. INenkova et al . 2008al ). 

There are two fundamental reasons for the different SEDs of clumpy torus models. First, 
the illuminated and back sides of individual clumps radiate at very different temperatures. 
This allows a certain location in the torus to have a large range of temperatures rather 
than the single temperature typical of continuous models. Second, and more important, 
gaps between the clumps allow a free, un-attenuated line of sight from the center to deep 
inside the torus. A clump at several r au b can emit at the local, distant dependent dust 
temperature spreading the high temperature regions over larger parts of the torus. There 
are two important consequences: the strength of the 9.7/rm absorption feature is much below 
its predicted strength in continuous distribution cases, and the anisotropy of the emitted 
radiation is smaller. An obvious limitation is the assumption of empty space between the 
clouds in a system where cloud collision and evaporation are unavoidable. 

The third type of composite models is, perhaps, the most realistic. In this case, the volume 
between clumps is filled with diluted dusty gas which absorbs part of the incident optical- 
UV radiation and part of the locally emitted NIR and MIR radiation. The main features 
of the pure clumpy case are preserved but additional attenuation by the inter-cloud dust 
and gas must be taken into account. F ig. [2] shows a schematic structure of the composite 
torus model of of Stalevski et al. 1 20121) and figure [3] exhibits several SEDs comput e d wit h 
this model. The SEDs are compared to the composite type-I SED of lMor fe Netzeii {2012). 

To illustrate the difficulty in deriving specific torus properties from comparing models 
with observations, I show in Figure [3] a very basic model that was constructed by combining 
the emitted fluxes of three (!) dusty clouds with the same column density, 10 23 cm -2 . The 
clouds are located at different distances from an accreting BH. Two of them contain pure 
graphite grains and the third graphite and silicate grains similar to what is found in the 
ISM (see details in the figure caption). This simplest possible theoretical SED is already in 
reasonable agreement with the observations suggesting that doubling or tripling the number 
of clouds is all that is required to get a very good fit. 

Finall y, we shou l d note two o ther types of phenomenological models, the tilted disk 
model of Lawrence fc Elvis (2010) and the combined torus plus dusty polar wind model of 


Honig ct al 


( 20131) . 


The first is discussed in 34] and the second in 
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Figure 3: 

Left: A comparison o f the lStalevski et al.1 (|2012h models (black lines) with the observed 
iMor fc Netz er (2012) composite type-I SED (red line) that contains dust emission from the NLR 
(see S- The models are marked with the viewing angle to the torus. Right: The SED of a three 
cloud obscurer. The column densities of the clouds are 10 23 cm -2 and the gas-to-dust ratio 
galactic. Two of the clouds contain pure graphite grains and are located at the graphite 
sublimation radius and three times further away. The third cloud contains ISM-type dust and is 
located at a distance which is 100 times further than the first cloud. The covering factors are 0.1. 
0.3 and 0.3, respectively. The solid black line is the combination of the three individual SEDs that 
are marked in dashed lines. The blue line is the SED of the central disk. 


2.3.2 Anisotropy and Covering Factor Each of the models discussed above has 
its own effective covering factor. For smooth torus models, this is simply a matter of geom¬ 
etry and the wavelength dependent optical depths. For clumpy models, it depends on th e 
number of clouds and their 3D distribution. In the straw person model of Antonucci (1993]), 
the covering factor is assumed to be the same for all sources of a given luminosity. The 
implication to large samples is translated to a simple expression: /c=N(type-II)/N(type-I). 
The realization of a significant spread in f c at a given lu minosity, suggests that the covering 
factor is a second unification parameter (lElitzurl 120121 ). In this case, the relative number 
of type-I and type-II sources still defines the mean fc but a source is more likely to be 
classified as type-II if its covering factor is larger than the mean and as type-I if it is smaller 
than the mean. Thus the high covering factor tail of the fc distribution contains more 
type-II AGNs, and the mean covering factor of type-Is is smaller than the mean of type-IIs 
if objects are selected by their optical properties. More quantitatively, if we define the 
covering factor probability distribution of sources with luminosity L as P(fc,L), and the 
probability for a type-I source from this population to have a covering factor fc, Pi(fc,L), 
we can write: 


P I ( fc , L ) = { l - f c ) P ( f c , L ) ; Pn ( fc , L ) = f c Ptfc , L ). (4) 
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Because of the finite probability of seeing the central source at any inclination through a 
clumpy medium (assuming no dusty inter-clump material), some sources classified spec¬ 
troscopically as type-I AGNs, have high inclination angles. Similarly, an obscuration by 
a single large cloud can result in the appearance of a type-II spectrum. This interesting 
possibility can be tested by observations. For example, the angle of polarization of hidden 
broad emission lines in face-on tori, is predicted to be v ery different from edg e-on tori. The 
presence of a lower density, dusty intercloud medium (IStalevski et al]|2012h changes this 
simple estimation but the overall dependence on the angular distribution of the clumps is 
still an important feature of all clumpy torus models. 

Anisotropy is a second characteristic of all torus models. This is most noticeable in 
continuous gas distribution tori where the smallest-distance dust, which is directly heated 
by the central source, can only be seen from small inclination angles. In clumpy tori, 
there are always gaps in the distribution and hot clumps can b e see n at larger incli nation 
angles. This reduces t he anisotropy. Inspection of several of the lStalevski et al.l (120121 1 and 
Nenkova et al. ( 2008bll models designed to fit observed SEDs, suggests that for wavelengths 


that are optically thick to the torus emitted radiation, the anisotropy is roughly proportional 
to cos i, where i is the inclination to the LOS. This is similar to the expected anisotropy of 
the central accretion disk. At longer wavelengths, where the dust optical depth is smaller, 
the emission is more isotropic and at FIR wavelengths the torus emits more like a sphere. 

A general, somewhat simplistic approach which does not take into account wavelength 
dependences is, nevertheless, useful for improving the estimated geometrical covering fac¬ 
tor for type-I AGNs. This approach requires an isotropy parameter b that varies from 1 
(complete isotropy) to 0 (maximum anisotropy), and a general dust anisotropy function, 
a(b,fc)- Assume a case where the (unknown) torus covering factor is fc and the total 
AGN bolometric luminosity, Lagn, is completely isotropic and seen, unattenuated, through 
the opening in the torus, i.e. (1-/c)Lagn escapes the system. The total dust emission is 
fc Lagn and the emission through the opening is Ld U st = (l — bfc)fc Lagn. Recall that flux 
measurements refer to luminosity per unit solid angle (i.e. divide the AGN and dust emis¬ 
sion through the opening by (l-/c)), we can define an observed ratio, R=Ld U st /Lagn and 
relate it to the covering factor by requiring R = a(b, fc)fc, where the anisotropy function 
is: 

1 “ bfc 

(5) 

This gives a(b, fc) = 1 for complete isotropy (6 = 1) and a(b, fc) = 1/(1 — fc) for maximum 
anisotropy (6 = 0), as required. For example, in a case where fc= 0.7, complete isotropy of 
dust emission gives R = fc, maximum anisotropy gives R = 3.33 fc, and an intermediate 
case with b = 0.5 gives R = 2.17 fc- 

Finally, in optically thick tori, the covering factor may also be related to the narrow line 
intensities since the smaller the covering factor is, the larger is the opening and more NLR 
gas is exposed to the central source radiation. As shown in § [5] below, this is a somewhat 
simplistic approach since the NLR gas does not necessarily fill the entire cone and since 
absorption by dusty inter-clump material close to the edges of the torus cannot be neglected. 
Moreover, the intensity of the narrow emission lines relative to the bolometric luminosity, 
Lagn, increase with decreasing source luminosity in a way which seems to be completely 
unrelated to the torus covering factor. 
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3 TORUS OBSERVATIONS 

The large number of successful IR space missions of the last 20 years, the on going suc¬ 
cessful X-ray missions, and the significant improvements in ground-based interferometry, 
revolutionized the study of AGN tori. They provide huge data sets that can be used to 
obtain reliable IR SEDs, allow multi-wavelengtli studies of many thousands of AGNs, and 
enabled to probe regions that are 1 pc in size near the BH. This section describes these de¬ 
velopments and their application to individual AGNs. The related redshift and luminosity 
dependences are discussed in m 


3.1 The Inner Torus Boundary: Dust Reverberation Mapping 

The variable luminosity of the central source results in time-dependent heating and ion¬ 
ization of the surrounding gas and dust that can be used to obtain valuable information 
about the geometry and physical conditions in these regions. Reverberation mapping (RM) 
has been used, for years, to map the d istribution of the line e mitting gas, mostly H/3, 


in the BLR (“line RM”, see Kaspi et al. 2000l . Bentz et al. 20131 . and references therein). 


Other experiments were designed to measure the time-dependent dust emission in the NIR, 


main l y in the K-band, in response to the variable continuum (“dust RM”, see e.g. Glasi 


1992, 


Suganuma et al. 20061. Kishimoto et al. 2007. Koshida et al. 2014. Pozo Nunez et al 


20141 1. The K-band wavelength, ~ 2.2/rm, is close to the peak emissivity of the hottest dust 


and the experiments provide information about the inner torus boundary 1 41.21) . 

The main results of dust RM experiments, now carried out successfully in about 20 low- 
to-medium lu minosit y AGNs, is a tight correlation of the form: — 0.4(L 46 ) 1/2 pc 

( Koshida et al I 20 H - where L 46 =Lagn/10 46 erg s 1 and I converted the measured V-band 
luminosity to Lagn assuming Lagn— 8ALa(5500)A. Thus, to a very good approximation, 
= (-Rsub.c) (eqn.[IJ. Note, however, that d ifferent value s for the dust sublimation 
radius have been used in the literature (e.g. Il-Coshida et al. 2014 ) which lead to suggestions 
that the measured V-K time-lag is 2-3 times smaller than expected from sublimation tem¬ 
perature considerations. The comparison of the H/3 and K-band RM results is shown in 
Figure [4] The measured K-band time lag is a factor 3-4 larger than the H/3 time lag, a 
result which is exp ected given what we know about gas and emissivity distributions in the 
BLR ( Netzer 20131 4 . 

RM experiments are not without difficulties. The dust transfer function (the dust emis¬ 
sion in response to a delta-function central source light-curve) must be complicated, and 
differs substantially from simple thin or thick shell. It depends on the torus and disk in¬ 
clination, the dimming of the disk rad iation towards its equato r ial plane, the occultation 
of some parts of the to rus, and more ( Kawaguchi fc Moril 2f)ld . Kawaguchi fc Mori 2011 


Pozo Nunez et al [20141 '). Nevertheless, the agreement between the dust RM results, and 


the predicted graphite sublimation radius, suggest that the inner torus wall is determined 
by the maximum grain temperature. 

Two comments about the relative sizes of the BLR and the inner torus walls are in order. 
The H/3 radius measured by RM experiments should not be mistaken for a single location 
BLR. In fact, the more successful BLR models req uire the ionized dust -free gas to extend 
all the way from the inner disk to the torus (e.g. Mor fe Netzer 2012h . In addition, the 
quoted distance ratio of 3-4 is only an average and actual ratios range from 1.5 to about 
10. This reflects the uncertainties on individual measurements and the likely large scatter 
in RM-based distances. 
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Figure 4: 

Various measures of the size of the central torus: MIDI-based measurements of R \/2 (red points 
and a red line), dust RM measurements based on V and K-band light curves and K-band 
interferometry (orange points and an orange line), graphite sublimation radius from eqn. |T] (black 
line), and the H j3 RM size (blue line ). The Typical 8 - 13//rn size is ~ 30 times larger than the 
typical dust RM size, (adopted from lBurtscher et alj d2013h . courtesy of Leonard Burtscher). 


3.2 Torus Properties: SED fitting 

The various torus models of Sj2] make specific predictions about the emitted SED which 
can be compared with NIR-MIR observations. This has been an active area of research 
because of several successful space missions, mostly the Spitzer space telescope ( Spitzer) 
and the wide-held infrared survey explorer ( WISE). The satellite data are supplemented 
with ground based, J, H, K, L and N-band observations with inferior wavelength coverage 
but much superior spatial resolution. The spectroscopic characteristics of the torus that 
allow to distinguish type-I from type-II SEDs are the 9.7 and 18pm silicate features and 
the l-8pm continuum slope. The comparison between observations and theory discussed 
below is biased towards phenomenological models since these are the ones that were used 
in almost all such studies. 

The spatial resolution of the IR observations is of great importance because of galactic 
scale emission that dilutes the characteristic signature of the torus. Stellar emission can 
contribute to the observed signal in the NIR and dust emission from SF regions can dilute 
the torus signature at 10-30 pm. In particular, strong polycyclic aromatic hydrocarbon 
(PAH) features, mostly at 6.2, 7.7 and 11.3pm, can contribute significantly to the MIR 
emission and dilute broad band photometric observations. In addit ion, dust emission from 


the NLR can contribute significantly to the 10- 30/im continuum IlSchweitzer et al 


2008 : 


Mor, Netzer fe Elitzur 20091 : Mor fc Netzer 2012 ). Because of this, small aperture spectro¬ 


scopic observations are preferred. Unfortunately, the Spitzer spectroscopic apertures are 
several arcsec in size and include significant c ontribution from the ho st galaxy in low lumi¬ 
nosity AGNs (see a large collection of those in Buchanan et al. ( 20061 ) ). A detailed compar¬ 
i son o f larg e and smal l apert ure spectroscopic observations is given in lAlonso-Herrero et, al. 
( 2011 ) and Lira et all ( 2013). _ 

The works of Ramos Almeida et al. 1 201 lh and Ramos Almeida et al. (2014) are good 
examples of high spatial resolution (~ 35 pc), very detailed observations of a very small 
number of objects. The results are similar to other works by this group and claim to be 
“clean” since they avoid galaxies with known dust lanes that affect the SED shape and 
the depth of the 9.7pm feature. The 1.5-19 pm SEDs obtained in those studies indicate 
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a considerably different slope for the two types of AGNs: For type-Is A Lx oc A 0 ' , and for 


type-I Is A L\ oc A 1 ' , with a rather large scatter in both slopes. However, the iLira et al. 


( 20131 1 work on a much larger sample of type-II sources, disagrees with some of the results. 
It shows a very large rang e of slopes, including some th at are similar and even flatter from 
the type-I slopes found bv IRamos Almeida et alj (120141 '). It seems that there is no canonical 
NIR-MIR slope for low luminosity type-II AGNs. 

Most attempts to derive torus parameters, such as the number of clumps along a certain 


direction and clump optical depths, give inconclusive results. Lira et ah (1201311 made a 


detailed co mparison between the rather different clumpy torus models of iNenkova et al 


( 2008bll and Stalevski et al. ( 2012l l and were able to obtain statistically satisfactory fits for 
both. One of the parameters of the fit is the inclination angle but, unfortunate ly , th ere is 
no correlation between the angles found by applying the two models. Lira et al. ( 2013h also 
found that SED fitting of about half of their type-II sources requires an additional very 
hot dust component, with dust temperature of ~ 2000 K, to explain the observed NIR flux. 
Such hot dust is not included in any of the clumpy models used in this work. It seems 
that the NIR-MIR spectra of AGNs of all types do not contain enough spectral features to 
distinguish one object from the next, and one clumpy torus model from another. In other 
words, torus models contain enough free parameters to fit almost any obs erved SE D. 
Th ere are several examples of SED fitting of luminous, type-I AGNs llSchweitzer et al. 


2008; Mor^^Jetze^^JShtzur 20091 : Roseboom et al. 20131 : Assef et al. 20131 : Lusso et al. 
2013 1. Schweitzer et al.l ( 2008h demonstrated that successful SED fitting of large aper¬ 


ture observations requires, on top of the standard torus model, two additional components: 
very hot, T> 1500 K emission f rom the torus and w arm dust emission from the NLR. 
This ap proach was later used by Mor fc Netzer ( 2012l l who combined the Nenkova et al. 
( 2008blf clumpy torus model with two additional components, one containing pure graphite 
grains and the other a dusty NLR. The hot dust component is not an integral part of the 


Nenko va ct al. (2008bj) models and its treatment is rather simple. After subtracting the host 


SF contribution, they fitted Spitzer spectra combined with NIR photometry and obtained 
satisfactory fits, like the ones shown in figure 0 for more than a hundred medium and high 
luminosity AGNs. The composite SED of this sample (red line in figure [3]) indicates a very 
flat, slope zero (in XL\) 2-20/j.m continuum and a turn down at longer wavelengths. The 
exact long wavelen gth dependence i s very sensitive to the assumed SF contribution. Indeed, 
a parallel effort bv lMullanev et all (1201 il l using a smaller sample of AGNs of both types, 
suggested a decline that starts at a longer wavelengths of about 30/rm. As shown earlier 
(figure [3|, a very simple combination of very few large clumps can fit such observations 
almost as well as the s ophis ticated clumpy models. 

The iMor fe Netzer ( 20121 1 SED fittin gs was used to put con straints on the various pa¬ 


rameters of the clumpy torus model of INenkova et al . (2008b). The results show a well 
defined inner radius, an ill defined outer radius, and a broad dist ributions o f cloud prop¬ 
erties. Many of these properties are similar to those found by ILira et al.l (120131 1 when 


fitting their type-II objects. The fittings also allow to compare the three major contribu¬ 
tions to the total dust emission. For the luminosity median of the population, L(hot torus 
dust):L(wa rm torus dust):L(NLR dust)«40:40:20. The relative NIR and MIR luminosi¬ 
ties in the Roseboom et al. ( 2013 1 sample, that does not consider the NLR contribution, 
are similar and somewhat luminosity dependent. The fits also provide an estimate of the 
luminosity dependent size of the NLR, 

n 0.67 


Rnlr — 460 


Tagn 


10 46 erg s“ 


pc. 


( 6 ) 
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Figure 5: 

Left: Observed (gray ) and SF subtracted (dot-dashed line) SEDs of two type-I AGNs from 
iMor &; Netzed ll2012l l. The diagrams show the three components used in the fit: hot pure graphite 
dust, a clumpy torus model with ISM-type dust, and dust emission from the NLR. The thick 
black curve is the sum of the three. Right: Composite type-I (q uasars) and type-II (radio 
galaxies) SEDs of z ~ 1.2 steep spectrum radio-loud AGNs from lHonig ct akl (1201 lf> illustrating 
the large spectral differences between the two groups. 


This dimension is more that two orders of magnitude larger than the silicate sublimation 
radius and is in agreement with the inner part of the NLR derived from direct imaging and 
and photoionization models. It corresponds to a mean grain temperature of ~ 160 K with 
peak emission at around 18/tm. 

An additiona l detai led study of the NIR-MIR SED of z ~ 1.2 luminous radio-loud AGNs, 
bv iHonig et al.l (1201 lfl . covers the rest-frame w avelength o f 1-17 /im. The sample includes 
3CRR sources that were originally described in Haas et al. 1 20081) and Leipski et al. 1 2Q10l) 
and are very similar in their IR-optical-UV properties to radio-quiet AGNs. Despite its 
small size (20 objects), the sample is most suitable for the comparison of type-I and type-II 
sources since it is chosen on the basis of the extended radio emission which is expected 
to be radiated isotropically. While the mean radio luminosity of the two groups is almost 
identical, the NIR-MIR spectrum is very different (see figure [5]). The type-I SEDs (r e ferred 
to as “quasars” in the diagram) are flat in \L\ and very similar to the Mor fc Netzer ( 2012 ) 
composite. The type-II SEDs (referred to as radio galaxies) are steeply rising from 1 to 7 
/im and are similar to some of the type-II SEDs discussed earlier. 

In summary, the typical NIR-MIR SEDs of type-I and type-II AGNs are rather different. 
Most, but not all type-Is show significant hot dust emission. Type-IIs show a much larger 
diversity with some objects similar to type-Is and others showing a steep decline from long 
to short wavelengths. Clumpy torus models provide adequate fits to many observed SEDs 
but seem to be degenerate with respect to their basic parameters. Because of this, models 
with different geometry, dust distribution and even inclination, can produce equally good 
fits to the same data. The inner region of the torus which contains pure graphite dust is 
clearly required by the observations but is missing from todays models. Given all this, it is 
not yet clear that current phenomenological models resemble, well, real AGN tori. 
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3.3 Anisotropy and Covering Factor 

As explained in §2.3.2, dust emission from optically thick tori is characterized by a wave¬ 
length dependent anisotropy. Several of the torus models described in fj2]are optically thick 
through most of the 1 20/rm wavelength range and SED differe nces between type- I and 
type-II sources can be used to derive the anisotropy. According to lHonig et al. (2011), this 
is the cause of the different SEDs shown in the right panel of figure [5] In this particular 
case, if the type-II inclination angle is 75 degrees, and the type-Is 39 degrees, the column 
density needed to produce this anisotropy corresponds to Ay ~ 150 mag. Even at such a 
large column, the optical depth is small and the emission is isotropic at A > 30/rm. 

A direct way to estimate the anisotropy is to measure the dust optical depth through 
sensitive search for broad Paschen and Brackett lines in the spectrum of type -II sources. 
The total number of objects studied in this way is small (see ILutz ct al. 120021 . and refer¬ 
ences therein) and the typical detection rate is only about 30%. The most useful results 
were obtained for type-II AGNs showing broad Balmer lines i n polarized l ight i ndicating 
a hidden BLR. The hand-full of broad Bra lines reported by ILutz et al . (2002) indicate 


that t(4.05 fj,m) ~ 3 mag., corresponding to column densities of order 10 23 cm~ 2 that are 
consistent with the columns derived from X-ray observations of the same sources. For sev¬ 
eral objects in this sample, with polarized broad emission lines but no direct broad line 
detection at 4.05/rm, the obscuration must be larger. As shown below in JfJ such a good 
agreement is not typical of many AGNs. 

The large NIR optical depths deduced from the Bra observations should not be confused 
with the observations of many type-Ii AGNs showing weak broad Ha lines and basically no 
broad H (5 indicating a few magnitudes of visual extinction. Such objects are common with 
numbers that are estimated to b e 10-20% of all low redshift AGNs (e.g. lLawrence fe Elvis 
19821 . Antonuccl 19931 . Elvisll2012 , see HTJ. They are thought to be reddened by galactic scale 
dust although in some cases, the source is probably seen through a small column density 
part of the torus. Reddening of high luminosity high r edshift AGNs is also a common 
phenomenon (e.g. Glikman et al. 2004 . Lacy et al . 20071 '). Such objects are found in IR 
and radio surveys and their weak, occasionally undetected optical emission suggests a large 
range of dust extinction, from that typical of many local type-Ii sources, to much higher. 

Their fraction in the general population must be high but is not accurately known 

As already alluded to, anisotropy influences covering factors estimates based on the mea¬ 
sured Ldust/L agn (eqn.(5j. The estimation of the mean and the distribution of fc requires 
an additional correction that takes into account the fact that a given Lagn is associated 
with a range of covering factors. In particular, a large fc increases the probability of 
classifying a source as type-II and the opposite is true for type-Is (eqn. |4j. 

Covering factors f or a large number of type-I AGNs have been obtained by using the 
Ldust /Lagn method (|Maiolmo_etjijJ|200I^ Trg^ei ^JCrolil^fc DuUemondl2008l : Mor fc Trakhtenbrot 


2011; Mor fc Netzerjl2012l; Roseboom et al.ll2013l; Lusso et al.ll2013h . The quality of the SED 

fits in such studies ranges from poor to very good, where spectroscopy-based fits provide 
more reliable measurements. Not all studies include the anisotropy correction term which 
results in overestimation of fc- All studies show a broad distribution of fc for a given 
luminosity, with a type-I population mean of 0.3-0.6. In general, there is clear indication 
for a decrease of the mean covering factor with increasing sour ce luminosity. Some of the 
studies attempted to separate hot (NIR) fro m warm (MIR) dust ( Mor fc Trakhtenbrotllioil 


Mor fc Netzer 20121 . Roseboom et al. 2013l b They show that the hot-dust covering factor 
varies over a larger range with some sources missing this component all together (“hot-dust 
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poor AGNs”, see Mor fc Trakhtenbrotl tl()\ il l). The fraction of such sources among lumi¬ 
nous AGNs is about 15%. In general fc (hot-dust) ~fc (warm dust) but the separation of 
the two is model dependent. A big disadvantage of this method as applied to IR-selected 
samples is the bias against low covering factor sources. 

The torus covering factor can be obtained dir ectly from modelin g the o bserved X-ray con¬ 
tinuum (see ©■ This has been demonstrated bv lUeda et alJ (120071 s ) and bv lBrightman fc Uedal 
( 20121 1 who made a careful comparison of the scattered and transmitted X-ray flux of a large 
number of type-II AGNs. They claimed to be able to measure scattered radiation which 
is as low as 1-2% and thus identify “buried” active BHs that are completely surrounded 
by geometrically thick tori with fc — 1- The work suggests that some 20% of high redshift 
sources in several deep X-ray fields are completely buried AGNs. 

An additional way to estimate fc is based on the observed intensity of the narrow Fe Ka 
line (SE3- Studies of large X-ray samples show that EW(Kq) decreases with increas¬ 


ing X-ray luminosit y, roughly like La gn 


but not with redshift (the “X-ray Bald¬ 


win effect” (see e.g. Chaudharv et al. 2Q10l . Ricci et~akl 20131 ')'). As explained, EW(Ka) 
depends on the combination of the torus column density, covering factor and metallicity 


(Ikeda. Awaki fc Terashimall2009l: Yaaoob et al.l2Qlol : Brightman fc Nandral2011a ; Ricci et al 


20131 '). The advantage of this method is that the selection of sources by their X-ray flux 
is independent of the measured quantity, EW(Ka). The possible complications are extra 
emission from Ka produced in broad line clouds and the fact that EW(Ka) depends also 
on the column density. Moreover, the measured LOS X-ra y column, and the mean column 
density of the torus, are not necessarily the same (see e.g. Brightman fc Ueda 2012 1. 

Finally, a complementary approach to find a population mean covering factor is to mea¬ 
sure the LOS obscuring column and associate small columns with type-I sources and large 
columns with type-II sources. The mean fc is then obtained from the ratio N(type- 
II)/N(type-I), where N stands for the number of objects. This method is discussed in 

m 

It is interestin g to examine the c overing factor of the most luminous AGNs. An example is 
the work of Weedman et al.l ( 2012 1 who compared the WISE (rest-frame 7.8/rm) and SDSS 
(rest-frame wavelength of 1350A) luminosities of the most MIR luminous type-I AGNs at 
redshift 1.5-5. The study assumed Ldust/L(7.8/rm)~ 3 and found extremely high values of 


Ldust, 10 4 


1 erg s , and fc> 1 for many sources. The problem with this and several 


other studies is the neglect of known selection effects and necessary correction factors, 
like the anisotropy correction term of eqn. [5] and the possible anisotropy of the accretion 
disk radiation. Moreover, these objects are probably unique among type-I AGNs in terms 
of their covering factor because they were selected to be the most IR luminous. Indeed 
the agreement between the luminosity distributions at 7.8/rm and at 1350A, is very poor. 
Finally, a large fraction of the sources may be reddened, because of the large dust covering 
factor, a possibility that was acknowledged by the authors. Taken together, it seems that 
the mean torus covering factor for the most luminous type-I AGNs is ~ 0.5. 

An important way to evaluate the combined effect of anisotropy and covering factor is 
to compare the X-ray and MIR luminosities, L(MIR)/L(2-10 keV). The comparison made 
at 6 and ~ 12/rm shows an almost identical mean ratio for the two AGN types, over four 
orders of magnitude in luminosity, with a typical scatter of about 0.3 dex ( Lutz et^L 2003 


Leverpon_et^lj200S l GanBbJ^OOojjAsmu^e^alJ^Olllj Honig et al.l2Qlol . Brightman fc Nandra 


2011bl . Gonzalez-Martfn et alJ 20131 . Merloni et al. 2014 ). This is a somewhat surprising re¬ 
sult given the covering factor and anisotropy discussed earlier. Naive expectations suggest 
that the larger covering factor of type-II sources will increase L(MIR)/L(2-10 keV) relative 
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to type-Is. On the other hand, anisotropy of the MIR emission would tend to decrease this 
ratio in type-II AGNs that are typically seen at higher inclination angles. Emission by NLR 
dust will affect large aperture observations and will contribute equally to the two types. 
This cannot be a significant source of confusion at 6/rm, where the NLR dust contribution 
is small, or when using small aperture observation that exclude most of the NLR emission. 
Is it possible that all these effects conspire to smear out the differences between type-I and 
type-II AGNs? There were several attempts to blame the unexpe cted constant ratio of 
L(MIR)/L(2-10 keV) on the method used to estimate L(2-10 keV) i Brightman fe Nandra 
2011bl . LaMassa et al. 2011 . Mavo fc Lawrencel 2013l f and to obtain absorption-corrected 
X-ray luminosities that take into account more complex clumpy tori, including partial cov¬ 
ering in some directions. However, it is difficult to see how all this can explain the above 
luminosity ratio.. 


3.4 Imaging Interferometry 

Perhaps the most important development regarding AGN unification is the significant im¬ 
provements in long baseline interferometry and the ability to resolve the central structure on 
a milli-arcsecond scale. The most significant results were obtained by three instruments: 1. 
MID-infrared interferometric Instrument (MIDI): The ESO interferometer working in the 
8-13rtm atmospheric window with projected baselines ranging from about 30 to 130 m (e.g. 
Jaffe et al. 2004 . Kishimoto et al. 2011bl . Burtscher et al. 20131 . Honig et al. 2013h . This 
instrument provided the largest sample. 2. The Keck K- band interferometer (KI) (see 
Swain et al. 20031 . Pott et al. 20ip| . Kishimoto et al. 2011a) with a projected baseline o f 
85m. 3. The ESO Astronomical Multi-BEam combineR (AMBER Weigelt et alJ 2012 1. 
Interferometric measurements allow to measure and model the size and shape of the IR 
emitting components on a very small spatial scale of less than 1 pc. The attainable reso¬ 
lution depends on the flux distribution in the source, the wavelength and baseline of the 
experiment, the redshift and the source coordinates (which determine the coverage of the 
uv-plane). The results are often quoted in units of -R 1 / 2 , which is the wavelength depen¬ 
dent radius containing half of the emitted flux. The MIDI latest results, summarized in 


Burtscher et alJ ( 20131 1 and illustrated in Figure [4j include 23 sources with 8-13/xm dimen¬ 
sions of 1-100 pc. The KI and AMBER samples are smaller but they provide valuable 
information about shorter wavelengths and hotter dust. 

The large scatter in the measured -R 1/2 at MIR wavelengths is very different from the 
very tight correlation between size and luminosity indicated by dust RM that probes regions 
closer to the inner torus boundary. Part of this must be related to the fact that in many 
cases, there are two distinct components in the resolved N-band interferometry. Translating 
this to a specific size is obviously model dependent and subjected to large uncertainties. 
This indicates a rather complex relationship between radius and dust temperature in the 
torus. Another surprise is the very similar size of type-I and type-II objects of the same 
luminosity (simple tori predict that type-II tori will look larger). Also, The measured K- 


band_I?i /2 is larger, by 0.2-0.4 dex, than the measured K-band RM size (IKishimoto et al. 


2011a| j. This is likely related to the different definitions of the two radii since R 1/2 is a flux 
weighted radius and Ri„ (or R S ub ) a response weighted radius. 

In a very small number of well observed nearby sources, the interferometry allows to 
obtain “model images” that reveal the rough structure, elongation, and inclination of the 
central dusty structure. As of 2014, there are four sources with this kind of information: 
NGC 1068 (type-II), Circinus (type-II), NGC424 (type-I) and NGC3783 (type-I). Of the 
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four, the spatial information obtained for the Circinus galaxy, a type-II AGN with a graphite 
sublimation radius of about 0.01 pc, i s perhaps the best. Th e emissivity maps of this source, 
in three MIR bands, were modeled bv lTristram et al.l (120141 ) and their combination is shown 
in Figure [6j The figure shows also the boundaries of the ionization cone based on HST 
imaging of the [O ill] A5007 line. The map shows two distinct components. The first is a 
thin (axis ration of 1:6 or less) ~ 1 pc disk-like structure, emitting about 20% of the MIR 
flux, with orientation which agrees with the orientation o f the smaller (0.1 — 0.4 pc) warped 
maser disk observed in this source ( Greenhill et al. 20031 1 . The second is a much larger and 
thicker elongated structure, roughly perpendicular to the first component. This elliptical 
shape component, elongated roughly in the polar direction, emits about 80% of the 8-13/rm 
flux with a clear wavelength gradient reflecting the dust temperatures in the various parts. 
None of this is in agreement with the shape and orientation of a simple central torus. 



Figure 6: 

A MIDI-based model-image of the Circinus galaxy showing the location of the maser disk 
(red-green-blue line), central thin MIR emitting structure (white), and a larger thick elongated 
structure with colors marking central wavelengths: red - 13//m, green - 10.5/vrn and blue - 8//rn. 
The thick red-green-blue elongated structure emits about 80% of the MIR radiation. Its 
orientation is very different from the expected orientation of a “standard” torus. The boundaries 
of th e ionization cone (on ly one side cone in this source) are shown in dashed blue lines (adopted 
from iTristram et al .1 J2014ll , courtesy of Konrad Tristram). 


A second exciting result is obtained by combining MIDI (N-band) and AMBER (K-band) 
measurements of the type-I AGN NGC 3783. The combined model-image indicates different 
elongations of the K and N-band isophotes. The most likely elongation of the K-band model- 
image (which is still somewhat uncertain) is in a direction parallel to the expected plane of 
the torus. The N-band model-image is of much larger dimensions and is elongated almost 
perpendicular to t he K-band model-i mage, along the expected polar axis. This geometry 
was interpreted by Honig et al. ( 2013h as a combination of an inner small torus, where the 
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dust is very hot, and a dusty polar outflow which, because of its distance, contains cooler 
dust. This is somewhat similar to the Circinus model mentioned earlier (where K-band 
information is not available). Interestingly, earlier ground-base d AO-assisted ob servations 
pointed out polar emission in a different source (NGC 1068, see Bock et al. 2000l) . 

While complicate d radiative transfer effe cts may result in the unusual model-images 
presented here (see Schartmann et al. 2005lh the interferometric results for Circinus and 
NGC 3783, and indications for similar structure in two other sources, suggest a possible 
change of paradigm. A central equatorial, small covering factor torus, may still be needed 
in order to explain the hot dust emission, but much of the cooler ~ 300 K dust may well 
reside in an additional component which has different geometry, size and optical depth 
and, perhaps, different dynamics (i,e, a wind component). A realistic model for this new 
structure must combine several different elements and take into account the fact that the 
hotter inner torus must also emit some of the observed MIR flux (e.g. 20% in the case of 
Circinus). It should also answer several outstanding questions about the clear pass to the 
ionization cone, indicated by strong unreddened narrow lines, if the flow is optically thick, 
or the lack of silicate emission in the Circinus spectrum if it is optically thin. One must also 
dismiss covering factor estimates that apply only to simple axisymmetric tori. In particu¬ 
lar, the assumption of /c(total)=/c(hot dust)+/c(warm dust) may not work in this type 
of geometry. Finally, it is interes ting to note that several hydrodynamics models suggest 
dusty polar outflows ( Wada 2012l l which may go part of the way towards in explaining the 
observations. 


4 ADDITIONAL AND ALTERNATIVE OBSCURERS 
4.1 Galactic-scale Obscuration 

Galaxy scale obscuration is common in many AGNs and numerous individual cases stud- 
ied by HST and ground based observa t ions have been documented s i nce the 1980s (e.g . 
Lawrence fe Elvi3 19821 ; Antonucci 19931 : Malkan. Goriian fc Tam 1998 : Prieto et al. 20141 ). 
Well known cases are some of the objects detected by optical spectroscopy and classified 
here as type-Ii AGNs. In several past publications they were referred to as narrow line X-ray 
galaxies (NLXGs). Some of the objects show kpc scale dust lanes, others might be affected 
by dust in the torus . The fraction of su ch objects might have been underestimated since 
many real type-IIs (IMerloni et all 120141 '). especially those with low quality spectroscopic 
observations, may belong to this group. As argued in §7, this group is a major source of 
uncertainty when coming to estimate the relative fractions of type-I and type-II AGNs 
Galactic obscuration is more common in low luminosity type-II AGNs. High spatial 
resolution observations of such objects reveal kpc scale dusty filaments that seem to be 
connected with much smaller dusty features close to the nucleus. This led iPrieto et al. 


( 20141 4 and others to propose that a central torus is not required for AGN unification at 
low luminosity. Large scale LOS absorption has also been found in several Compton thick 
AGNs studied with Spitzer (jCoulding et al. 2012j). In this sample, AGNs in face-on galaxies 
show little 9.7/rm silicate absorption, as predicted by clumpy torus models, and AGNs in 
nearly edge-on galaxies, or in mergers, show much stronger silicate absorption indicating 
dust in the galaxy. Thi s adds to previou s similar results obtained for Compton thin type-II 
AGNs 1 Deo et al. 2007 . Deo et al. 20091 4. In some cases, galactic ex tinction is so large that 
it can completely obscure the central AGN at optical wavelengths I Colliding fc Alexander 
2009l l. The connection, if any, to the heavily reddened AGNs mentioned earlier, is not yet 
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clear. 

Supplementing the detailed investiga tion of individual sources, are lower spatial resolution 
studies of many thousands of AGNs like L agos ct al . (201 11). This study finds that spiral and 
elliptical type-I hosts are strongly skewed towards face-on galaxies, while type-II hosts are 
skewed towards edge-on orientations. The paper also shows that Galactic dust absorption by 
itself cannot explain the absence of broad emission lines in type-II sources. The inclination 
of the host galaxy strengthens the idea that the angular momentum of the material which 
feeds the black hole retains a memory of its origin in the galaxy, i.e., statistically speaking, 
galaxy disks, central tori and inner accretion disks, are generally alig ned. This issue, that 
was investigated in numerous papers, is still unresolved. For example, Kinney et al. ( 2000l ) 
used information from radio jet position angles and claimed to find no correlation between 
jet direction (i.e. inner accretion disk orientation) and the disk of the host galaxy. 


4.2 X-ray Obscuration by Gas in the BLR 

The simple idea of a single torus-like obscurer suffers from various limitations, especially 
when it comes to very large column densities. For example, the observed optical red¬ 
dening is much smaller than the value obtained by translating the measured X-ray col¬ 


umn to foreground extinction assuming galactic dust -to-gas ratio (e.g. iMaiolino et al 


2001 


Risaliti, Elvis fe Nicastrd 20021 : Goulding et al. 20121 ) . There are also X-ray variations that 
indicate large changes in the obscuring column on very short time scales (a few days or less), 
which is not consistent with the idea that the obscuring material is outside the sublimation 
radius. This raises the possibility that obscuration by dust-free gas inside the BLR affects 
the X-ray absorption and hence also the AGN classification. 

The typical dimensions of BLR clouds are estimated to be 10 13-14 cm and their typical 
column density is of order io 23 ~ 24 cm -2 . This is roughly the size of the central X-ray source 
in low-to-intermediate luminosity AGNs (a few R g ). For an X-ray continuum source ~ 10 14 
cm across, and a cloud velocity of order 3000 kms' 1 , the occultation time is of order 4 
days. Can such clouds occult the central source giving the impression of a variable X-ray 
source? Recent observations and modeling suggest that this is indeed the case in almost 
every object where the observations are of long enough duration to discover such “eclipses”. 

There is now evidence from very detailed studies of a small number of type -I and type-II 
AGNs, for eclipses by large column density blobs ( Risaliti et al. 20071 . Eoi lP). In the best 
studied case, NGC 1365, the X-ray sour ce underwent a tot al eclipse in less than 48 hours 
and then emerged within the same time (IRisaliti et al.ll2007l ). Detailed modeling of this and 
similar events, by Maiolino et al. ( 20101 ) ■ provided information about the likely location, 
velocity and size of the occulting gas, a ll cons istent with being part of the BLR and inside 
the sublimation radius. Maiolino et al. ( 2010 ) speculated that the occulting clouds might 
have developed cometary tails, which are manifested as a gradual decrease in the column 
density as clouds emerged from the eclipse. Different ti me dependent X - ray a bsorption 
events from inside the sublimation radius are described in Turner fc Miller ( 20091 ). 

There are more studies that are no t as detailed but include a larger num ber of AGNs 
and occultation events. The work of iMarkowitz. Krumpc fc Nikuttal (J2014) is based on 
long term monitoring of about a dozen AGNs of both types. These objects were monitored 
continuously by the Rossi X-rat timing explorer (RXTE) and the total duration is of order 
10 years. Given the known BH mass in these sources, one can work out the expected crossing 
time of dusty and dust-free blobs of certain properties at different locations. The results 
suggest that many of the clumps are, indeed, inside the sublimation radius with properties 
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that are consistent with being BLR clouds. Other blobs are larger and move slower with 
locations that are consistent with being inside the dusty torus. Interestingly, in several of 
the cases the clump location, ass uming Keplerian orbits, seems to be clo se to the maximum 
allowed dimension of such tori. Markowitz. Krumpc fe Nikuttal (j2014j) presented also an 
analysis of the non-variable X-ray absorbing columns and claim they are in conflict with 
the predicted dimensio ns, column density and number of clumps in clumpy torus models. 
An additional study, by Torricelli-Ciamponi et al. ( 20141 1 , includes more sources but shorter 
duration monitoring sensitive to events that last a few days at most. It presents numerous 
events associated with changing in hardness ratio of the X-ray continuum, all consistent 
with being caused by transiting blobs inside the BLR. The observations are not detailed 
enough to allow a robust determination of the column density of the absorbers. 

The X-ray observations provide a beautiful illustration of the connection between dust- 
free and dusty clumps in one general region that contains the BLR, closer in, and the 
torus, further away from the central BH. The earlier apparent disagreement between optical 
reddening and X-ray column density disappears when taking into account the combination 
of inner dust-free and outer dusty absorbers. The observations also help to explain the 
nature of some of the largest obscuring columns found in AGNs. 


4.3 Tilted Disks 


Is an axisymmetric torus-like structure the only type of a central obscurer that can explain 
the various aspects of AGN unification? Can other structures help to resolve several of 
the known disagreements between simple torus models and observation s? An interest ing 
scena rio is based on the idea that a central disk can replace the torus (ILawrence fe Elvis 
2010). This ad-hoc phenomenological model is based on the idea that the angular mo¬ 
mentum of the external inflowing gas is initially mis-aligned with the axis of the central 
accretion disk. The mis-alignment is kept down to very small distances, of order one pc, 
and the inflow takes the form of an in clined, large dusty geometrically thin disk, which 
may also be warped (see Caproni et al. ( 20061 1 for general mechanisms resulting in warped 
disks). The large disk presents a wide solid angle to the inner disk and part of the UV- 
optical radiation i s abso rbed and re-emitted in the IR. This simple geometry was used by 
Lawrence_& - j31vi^J 201C ) to predict the covering factor distribution which was later studied 
by Roseboom et alJ ( 201 3) using WISE data of a large type-I sample. There seem to be 
a reasonable agreement between model and observations but the predicted is similar to a 
normal distribution, which is expected also in other scenarios. 

The main virtue of the tilted disk idea is the assumed mis-alignment between different 
accretion episodes onto the BH. This is expected in accretion events whose final stage take 
place in a region which is 2-3 orders of magnitude smaller than the region where the flow 
originates. Unfortunately, the exact geometry has not been shown to be consistent with 
realistic hydrodynamical simulations, the location of the sublimation surface has not been 
calculated and compared with RM results, and there is no comparison with IR interferom¬ 
etry, IR SED and high quality polarization measurements. 


5 IONIZATION CONES 


An axisymmetric central obscurer with a central opening results in an anis otropic two- 
sided conical shaped ionizing radiation pattern (“bicone”). Early obs ervations (IPoggeill988l : 
Tadhunter fe Tsvetanov 19891 ; Mulchaev. Wilson fe Tsvetanov 1996h identified such cones 
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in a number of nearby systems and used them to deduce the torus inclination angles. The 
most detailed early studies are those of the archetype source, NGC 1068, and the most 
detailed mapping utilized the strong [Oiii]A 5007 line (the “[O ill] bicone”). The more 
recent studies include larger samples, more detailed imagining in the light of more emission 
lines, and additional kinematic information. The mapped emission lines now include low 
ionization species like [Ol], intermediate ionization ions like [O III] and [Oiv], and very 
high ionization ions like [Sivi]. High resolution images obtained with the Chandra X-ray 
observatory ( Chandra) allow mapping such bicones in individual X-ray lines like [O vn]. 
The most detailed observations of today are based on HST imaging, HST lo ng slit 


spectroscopy, and ground based integral field unit (IFU) studies ( e.g. Fischer et al. 120131 : 
Muller-Sanchez et ahlEoil : Riffel. Storchi-Bergmann fc Wineelboi.ll . and references therein). 


The number of sources mapped in detail is more than 50 and those with reliable kinematic 
information more than 20. The highest redshift is about 0.05. The observations are detailed 
enough to allow photoionization modeling of the gas in different parts of the NLR includin g, 
in rare cases, following the metallicity gradient across the source (e.g. Dopita ct al.ll2014l h 
The ability to distinguish low and high ionization regions across a spatially resolved NLR, 
provides evidence for ionization outside of the nominal [O ill] bicone. In several well studied 
cases, there is evidence of “filtering” of the ionizing radiation in the parts producing the 
strongest low ionization lines by clumpy, ionized absorbers located tens of parsecs or less 
from the nucleus, consistent with the torus dimension. A nice visual evidence is provided in 
NGC 4151, where HST imaging shows that the [O ml A50 0 7/1Q ill A3727 line ratio is lower 


near the edges of the bicone compared with its axis (IKraemer. Schmitt fe Crenshaw 


Such filtering was suspec ted i n earlier works bu t without the information about the location 

fconnh . 


of the filtering material (Alexander et al 


Line radiation from ionization cones provide ways to improve the estimate of the bolomet- 
ric luminosity of type-II AGNs where the non-stellar continuum is not directly observed. In 
particular, the [O ill] A5007 line was found to be a good l uminosity indicator b ecause the line 
is strong and not affected much by stellar absorption (IHeckman fc Bestl 2014j). This method 
is not without difficult ies. T he line is not a very good luminosity indicator in low ionization 
AGNs like LIN ERs dNetzeill2009ll and it is now known that L(fO ill] A5007) /Lagn is chang¬ 
ing with Lagn ( Stern fc Laorl 20121) . Line reddening is an additional uncertainty.. In par¬ 
ticular, reddening by dust outside the NLR depends on the source geometry and torus incli¬ 
nation to the LOS and is, therefore, different for type-I and type-II sources (line attenuation 
by dust inside the NLR clouds does not depend on these factors). Indeed, observations show 
that the emission pattern of some of the lines is more isotropic than in others which leads to 
the suggestion t hat ]Q iv] A25.9)im is the best luminosity i ndicator because MIR reddening 
is much lower (I Kraeme r et al. 201 ll : IWeaver et al. I I2OI0I: iDiamond-Stanic. Rieke fe Rigb-d 
20091 : Dicken et alJ 2014h . In fact. I([0 ill] A5007)/I([0 iv] A25.9pm) is typically lower by a 
factor of ~ 2 in type-II sources probably due to inclination-dependent reddening. Such a 
situation would naturally occur if the dust causing the extinction is just outside the [O ill] 
bicone and the bicone axis is closer to the plane of the sky. 

Soft X-ray mapping of several NLRs supplemented by Chandra and XMM-Newton high 
resolution spectroscopy, clearly indic ate X-ray line em i ssion, mostly Ovil570e!7, that gen¬ 
erally overlaps with the [O ml bicone (Sako et al.ll200(ll: Ogle et al . 200f3 : iKinkhabwala et al, 

i3 


2002 : Bianchi. Guainazzi fc Ghiabergel 20061 : IWane et al.ll2011ah . A s pectacular e x ample , 


showing X-ray emission up to about 7 kpc, in Mrk573, is given in IPaggi et alJ (20121). 


The LOS attenuation of X-ray lines depends only on column density and metallicity, 
and not on the dust content in the gas. Therefore, observations of X-ray lines pro- 


28 Netzer 






































































vide an additional method to measure this column. Kraemer et al. ( 2011 1 analyzed the 
O vn 570 eV/[0 iv] A25.9pm line ratio in a sample of about 15 sources and claimed an at¬ 
tenuation of the X-ray line by a column density which is consistent with the derived amount 
of reddening assuming ISM dust/gas ratio. While not a direct evidence, this is consistent 
with the idea of a dusty ionized gas outside and further away from the [O ill] bicone. Inter¬ 
estingly, X-ray emission is detected also in regions perpendicular to the optical bicone. This 
was interpreted as leakage through a clumpy, r elatively low column density absorber, which 
may be the torus. iKreimever fe Veillemd ( 20131 1 used optical line ratios, to compare the flux 
emerging from the opening of the cone to the flux leaking through the torus. According to 
this work, the leakage amounts to 30-50% of the radiation emitted by the central source. 
This fraction is very high and confirmation in other sources is required. 

The kinematic stu dies provide information about the motion of t he line emitting gas 
inside the cones (e.g. Muller-Sanchez et al. 20 111 . Fischer et al. 2014 1. Fischer et al. (2013) 
presented multi-slit observations of 17 AGNs, 12 of which are type-IIs, with clear signatures 
of biconical outflows. The velocity maps allow fits to simple kinematic models which are 
detailed enough to infer the inclination of the NLR to the LOS. In general, the angle 
of inclination is larger in type-II AGNs consistent with the view that these sources are 
seen edge on. Interestingly, there is no correlation in this small sample between the NLR 
inclination angle and the orientation of the disks of the host galaxy (see both supporting 
and contradicting results in earlier sections; t his issue is not yet se ttled). 

Kinematical modeling of ionization cones by Fischer et al. ( 2013) provides new invaluable 
information about ionization cones. It reveals a unique pattern of gas outflow with rocr 
out to a certain radius, about pc, and slowing down to zero beyond this radius. 

The fitted models are good enough to constrain the bicone inclination to the LOS and, 
in some cases, its opening angle. Another interesting finding is a strong increase of X-ray 
column density with increasing angular distance from the bicone axis, regardless of whether 
the column is neutral (type-IIs) or ionized (type-Is). For the type-I sources, where the 
X-ray absorbing gas is highly ionized (a warm absorber) it suggests a connection between 
the inner walls of the torus and the origin of the X-ray ionized material. The study also 
suggests a correlation between mid-IR color and polar angle (see For a handful of type-I 
sources, there is a hint for increasing FWHM(H/1) with polar angle. If verified in larger 
samples, it would indicate a rotational component for the BLR gas. Unfortunately, the 
results cannot be applied to the entire population partly because only one in three sources 
had good enough data to be modeled in this way, and partly because in some of the sources 
the m odel is based on a single slit observation. One good case is shown in Fig. [7| _ 


ThelMiiller-Sanchez et al.l (1201 il l work supports some but not all the findings in 


Fischer et al 


( 2013l l. This work provides imaging and kinematics of the coronal line region, the innermost 
part of the NLR where the strongest lines are of very high ionization, e.g. [Sivi] A1.96pm 
with ionization potential of 167 eV. In this region, the gas is probably dust-free. The study 
of seven sources shows a complex dynamics of the gas including a rotating disk and an 
outflow bicone. There are agreements and disagreements between the two works. This is 
illustrated in the case of NGC 3783, a well known type-I source with a large column of LOS 


ionized gas and interferometric data suggesting dusty polar outflow (see SJ3J. iFischer et al. 


( 20131 1 modeled this source using HST imaging and a single slit position. They found 
s mall inclination of 15 degrees to the LOS, as in a type-I AGN. The IFU measurements 
of Miiller-Sanchez et all ( 2011 1 resulted in a kinematic map based on [Si vi] A1.96pm and 
other coronal lines in a region very close to the base of the ionization cone. It suggests a 
large inclination, close to edge-on view, with the cone opening cutting into the plane of the 
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Figure 7: 

Kinematic mapping of the type-II source Mrk 34 from iFischer ct al.l (120131 ) showing slit positions 
(top left), velocity, flux and FWHM of [O III] A5007 in different positions (right panels), and the 
resulted bicone geometry (bottom left). The green and and blue points represent data obtained 
with the G430L and G430M gratings, respectively. 


host galaxy. If correct it means that the broad emission lines are seen through the (leaky?) 
torus. 

In conclusion, the study of ionization cones provides valuable information about the 
central obscurer in AGNs. The detailed mapping in the light of numerous lines confirms the 
general bicone shape and shows well defined borders between ionized and neutral gas. It also 
reveals the kinematics of the gas and the overall inclination to the LOS. The inferred models 
are more concrete and better justified experimentally than the rather vague information 
about the torus provided by fitting NIR-MIR SEDs. 


6 THE TORUS GALAXY CONNECTION 


General scenarios linking BH feeding and SF a ctivity have been descri bed in numerous 
papers and reviewed recently in this journal by [Heckman fe Best! (120141) . These concern 
global galaxy properties and apply to SF activity in the central kpc and beyond. They do 
not address the more direct connection between the two types of activity very much closer 
to the central BH. What is known about the central 10-100 pc of nearby AGNs? What 
is the stellar population in this region, the atomic and molecular gas content and the SF 
activity? Is there evidence for a large central disk and is it connected to the torus? 

Detailed observations of the central regions in nearby AGNs s hed light on s o me of these 
questions. Adaptive optics (AO) assisted observations, like in I Davies et al . (120071 ). can 
probe regions that are 10-30 pc in size. As of 2014, th ere are about 15 sources observed 
at this resolution. Another work, by Esauei et al. ( 2014! ). uses seeing limited N-band spec¬ 
troscopy in 29 sources and probes a region of ~ 65 pc in size. The aim is to investigate 
SF near the nucleus using the 11.3/rm PAH emission feature. About 45% of the objects 
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in the sample show this feature and stacking the spectra of the remaining sources suggests 
its presence in many others. There are several other works usi ng seein g lim ited optica l 
observations, with a typical spatial resolution of 100-300 pc (e.g. Riffel et al. 2010| . 2011 '!. 
Important questions related to the conditions in a region which is slightly larger than the 
size of the torus, or the sphere of influence of the central BH, are related to the age of 
the stellar population and the SF rate (SFR). As explained below, there is no consensus 
regarding these issues. 

Standard SF indicators are PAH features, optical absorption features and hydrogen re¬ 
combination lines. The recombination lines are the best direct indica tors for on-go i ng SF 
(wh ich is defined here a s the presence of O-type stars). According to iDavies et al.l (120071 1 
and lDavies et al.l (120141 1. the equivalent widths of the Br~, lines in the 15 AGNs studied by 
this g roup, practica l ly rul e out on-going SF and indicate post SF population. On the other 
hand. lEsauei et al.l (12014 1 claim to detect SF activity based on the detection of PAH emis¬ 
sion features. This apparent contradiction can be understood by noting that PAH emission 
is not necessarily related to the presence of O-type stars. In fact, there are indications that 
PAH emission is more sensitive t o B stars than O stars, i.e. to an environment that forme d 
stars in the past 100 Myr or so ( Peeters. Spoon fc Tielens 2004 j_ DIaz i Santos_et_aI. 2Q10l l. 


Moreover, the calibration of L(PAH) versus L(SFR) us ed by 


the SFR is based on work by Diamond-Stanic fc Riekel ( 2012 1 and Rieke et al 


Esguej_e£jji (2014 ) to d erive 


2009ll who 


used templates of SF galaxies. This uses SFRs and PAH luminosities obtained through 
large apertures that occasionally include the entire galaxy, and is equivalent to integrating 
over ~ 100 Myr of SF history. 

Returning to AGNs, the observed PAH emission close to the center is entirely consistent 
with the idea that the presently observed BH accretion episode is a few Myr old and BH 
accretion might have started after all the O-type stars in the region ended their life. In 
fact, the onset of AGN activity might have terminated all SF in the immediate vicinity 
of the BH which is usually categorized as “AGN fee dback” ( Davies. Mark fc Sternberg 
20121 : Krios et al. 20 111 : Hailev-Dunsheath et al. 2012l l. Additional support for the post 


SF environment, and occasionally an even older stellar population very close to the centra l 


BH, is obtained from oth er stu dies of the 100-300 pc central region ( Riffel et al. 2Q10l . 2011 
Storchi-Bergmann et al. 20121 . Schnorr-Muller et al. 2014). 


Some observations have sufficient spatial and spectral resolution to derive the mass an d 
age of the stars in the inner 50 pc. For the few cases this was done (e.g. IDavies et al . 2007]), 
the velocity field is dominated by random motion typical of a spherical system and the 
stellar mass is of order 10 s Mq. This mass is an order of magnitude larger than the mass of 
the central BH in questi on. These observations set the stage for the stellar-wind model of 


Schartmann et al 


q uesti c 

120101) 


discussed in m The idea of BH feeding by stella r wind material 
has also been explored in connection with LINERs ( Heckman fc Best! 2014, and references 
therein). The argument in this case is not based on spatially resolved observations but 
rather on the extremely low Lagn / L^dd in LINERs. It is not at all clear whether the two 
independent scenarios have anything in common. 

Studies of molecular gas in the vicinity of the torus do not yet have the required spatial 
resolution to resolve the torus, or even its close surrounding, and hence were given low 
priority in this review. The co ming ALMA observations will, no doubt, change this situation 
(see review by Maiolino 2008 :1. It is already clear that high concentrations of molecular gas 
on ~ 100 pc scales are seen in many AGNs with column densities that are large enough to 
confuse the pure torus absorption picture. There are clear indications for outflows that are 
possibly associated with the direction of the bicone, but no indication of inflow. Some of 
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these studi es are detailed in Hickset_aL ( 200flh . Oarcfa-Burillo et al. (2014), Combes et al, 

( 20131 4 and Aalto. S. et al. (l2014ll where references to older observations can be found. 

None of the above observations and theoretical considerations can be applied, by simple 
scaling, to higher BH mass, more luminous AGNs. First, present day observations do 
not have the spatial resolution needed to investigate regions of this size in such objects. 
Second, feeding 10 8-9 Mq BHs in high ionization AGNs by stellar winds would require a 
much larger stellar cluster and may involve very different conditions. Direct mass inflow 
from the host is l ikel y to play the central role in such cases. Finally we note recent results 
from I Merloni et al.1 (120141 4 that compared BH accretion rate and torus properties with total 
stellar mass and SFR of the host galaxy. They indicate that the kpc scale stars and gas 
do not seem to have a direct causal connection with th e BH and its im mediate vicinity. A 

1 201lh . 


similar lack of correlation at low redshift is shown in Koss et al 


7 UNIFICATION OR EVOLUTION? 

7.1 Luminosity and Redshift Dependencies 

As explained, there are two categories of methods used to estimate fc■ The first category 
include all methods for measuring the covering factor of individual objects (see ©■ They 
include the Ldust/LAGN method, the EW(Ka) method, and the method based on the scat¬ 
tered X-ray continuum. The Ld us t/LAGN method is the one resulting in the largest number 
of measured fc- The second category includes methods that apply to the mean covering 
factor of the population. This is obtained by classifying a source as either type-I or type-II 
AGN, and using the resulting N(type-II,L)/N(type-I,L) as an estimate of the mean value 
of fc- One such method is based on the measured X-ray absorbing column and the AGN 
type is determined relative to a fiducial absorbing column that separates the groups, (e.g . 


Uedu i- e^ - alj2003 ]_ Gilli. Gomastri fc Hasingeill2007 ; Hasingerll2008l : Brightman fc T Jedall20 1 ‘A 

Malizia et al.ll2012l '). Another is based on counting the number of broad and narrow emission 

line AGNs in optical and IR selected samples. Ideal l y, all methods should agree but this 


is certainly not the case (e.g. Lawrence fc Elvis 201Cll . Mayo fc Lawrencel 2013 . Lusso et al 


20131 . 20141 11. The reasons range from the poorly understood physics of the torus, to vari¬ 
ous selection effects, uncertainties in estimating Lagn, intrinsic reddening, and more. For 
example, the La us t/LAGN method is based on dust emission from the torus and has little 
to do with BLR properties. However, X-ray absorption by gas in the BLR can increase the 
total measured X-ray column thus decreasing Ld us t since the ionizing radiation absorbed 
by the dust-free gas never interacts with the dusty torus. Here I focus on the large sample 
approach, in particular on X-ray based studies. 

The hard X-ray (E> 2 keV) spectrum of most AGNs can be described by a simple power- 
law, N p h oc E ~ t , where N p h is the number of emitted photons per unit energy and time. 
A remarkable feature is the very small scatter in the power-law slope, F 2 _iokeV = 1.8 ±0.2. 
This allows us to estimate the obscuring column density even in poor X-ray observations 
provided we are looking through a single absorber. All that is required is a measure of the 
hardness ratio, e.g. F(0.5-2 keV)/F(2-10 keV), and the redshift. The resulting accuracy in 
measured column is quite remarkable (~ 0.3 dex), with even higher accuracy in cases of high 
photon flux that allow better modeling of the intrinsic continuum. The hydrogen column 
densities derived in this way range from very thin (about 10 21 cm -2 ) to Compton thick 
(> 1.25 x 10 24 cm -2 ). Some of the results obtained so far are hampered by selection effects. 
In particular, type-II sources are more likely to drop below the flux limit of the samples 
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because of the increased absorption. T his is less noticeable , but still present, in samples 
selected by their very high energy flux ( Malizia et al. 2012i l. The hardness ratio method 
can fail, sometimes dramatically, if the absorption is due to two or more LOS components, 
such as a combination of BLR clouds and a d usty torus gas (§4.2) 

A very detailed study of X-ray obscuration is lMerloni et al.l (120141 1. This work is based on 
observations of the two-degree XMM-COSMOS field and includes 1310 sources of all types 
selected on the basis of their absorption corrected 2-10 keV luminosity in the redshift range 
0.3-3.5. The rest-frame 2-10 keV flux limit is 2 x 10~ 15 erg s -1 cm~ 2 . The study supersedes 
earlier work of this type (e.g. Ueda et al. 20031 : Gilli. Comastri fc Hasinger 2007 : Hasingei 
2008.; Malizia et al. 2012h that used data from other fields and in struments. I t suppl ements 


a comprehensive study of type-I sources in the s ame field by lLussoet_alJj 2013 ). The 


spectr al properties of the sources are described in Brusa et al. (|2Q10l l and iMainieri et al 


(20LlJ). The complete redshift information is based on spectroscopy and photometry and 
is of crucial importance since it allows to use a simple K-correction and thus circumvent 
the well known bias of flux limited samples that tend to pick up more obscured sources 
at higher redshifts. A potential complication is a large number of Compton thick sources 
that are unaccounted for. This number is still uncertain and depends on the survey, energy 
range us ed, and the spectro scopic analysis. It is generally e stimated to be ~ 20% at low 
redshift ( Burlon et al. 12011 ) and ~ 40% at high redshift (e.g. Brightman fe Ueda 2012l b In 
the Cosmos sample which is not very deep, the numbers are smaller. There are also issues 
regarding the best way to correct the X-ray luminosity for absorption, in particular the 
simplistic way based only or hardness ratio with its well known limitations. 

An independent method to verify the estimated numbers of obscured and unobscured 


AGNs is based on mode ling the spectrum of the cosmic X-ray background ('e.g. lComastri et al 


1995, Ueda et al. 20031'! . I will not go i nto the details of this method w hich was reviewed 


and discussed in numerous papers (e.e. ICilli. Comastri fe Hasingerll2007lL 

An important question is whether type-I and type-II sources classified by their optical 
properties can be distinguished, unambiguously, by their X-ray obscuring column? The 


answe r is definitely no although type-II are, on the average, more obscured. Merloni et al 


( 2014 ) introduced a dividing column density of Nh — 10 21 ’ 5 cm -2 to separate the two 
types. They found about 30% optically classified type-I sources with obscuring columns 
larger than this value and about 30% of optically classified type-IIs with obscuring columns 
smaller than this column. In fact, many, perhaps most of the objects in these two groups 
may have little to do with obscuration by the torus. The sources in the latter group are the 
real type-II AGNs discussed in earlier sections, and at least some of the objects in the first 
group are type Is with an X-ray source occulted by gas in the BLR or outflowing ionized 
gas in the opening of the torus.. These two sub-groups were not fully appreciated by the 
old unification scheme and are definitely influencing the N(type-I)/N(type-II) statistics. 

Many of t he real type-II AGNs in the local universe are of v ery low luminosity and 
I/AGN/I/Edd ( Brightman fc Nandra 2011bl . Marinucci et all 2012l b Those in the XMM- 
COSMOS sample, at redshift of ~ 0.3, are more luminous and can be considered inter¬ 
mediate luminosity AGNs. An important clue to the nature of these sources is the presence 
or absence of a central obscurer. This can be tested by exami ning the L(MIR )/L(2-10 keV) 
relationship discussed in (J3] The real type-II objects in the lMerloni et al.l (2014) sample 
do not differ in this respect from other type-IIs. There is also direct evidence from Spitzer 
spectroscopy that the MIR SEDs of these sources are simi lar (b ut not identical) t o the SEDs 
of type-II sources showing hidden broad emission lines ( Tommasin et al. 2O10l l. It seems 
that many real type-II AGNs do not have enough high density gas to result in a detectable 
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BLR but have central tori. This raises an interesting question about the relationship be¬ 
tween the central source luminosity, the inner dimension of the torus, and the presence or 
absence of broad emission line gas. This question is discussed in m below. 

Returning to the type-I sources with X-ray absorption, we note that for a Galactic dust- 
to-gas ratio, the above dividing column corresponds to Ay ~ 1.7 mag. If the obscurer is 
part of the dusty torus, then the intrinsic optical-UV spectrum of 30% of all type-I objects 
are expected to be heavily reddened. This is definitely not the case. As noted earlier, the 
obscuring gas may be part of the BLR in which case it is dust-free. Alternatively, the 
absorbing gas is highly ionized outflowing, dust-free material (“warm absorber” which is 
known to be common in AGNs ( Turner fc Milleill2009l b The location of this gas maybe well 
outside the BLR, including in the opening of the torus. 

Figure [8] shows the fraction of obscured AGNs defined in two different ways as a function 
of L(2-10 keV). The dependence on X-ray luminosity of the two groups is very different, 
showing a steep decrease in obscuring fraction if classified by optical properties, and a 
very flat and weak dependence if X-ray classification is used. Clearly, the large number of 
real type-II AGNs is the main reason for the difference between the two methods at low 
luminosity. This suggests a new, perhaps more intuitive definition of AGN types based only 
on the presence or absence of a central LOS obscurer. According to this classification, real 
type-IIs belong in the class of type-Is (unobscured AGNs). In this case, the relative fraction 
of the two groups up to redshift of about 1.5, is about 1:1 according to both the X-ray and 
optical methods. 



Figure 8: 


The fraction of obsc ured and unobscured AGNs, at various redshifts and L(2—10 keV), from 
iMerloni et alj ll2014l '). The left panel shows the fractions obtained by using optical classification. 
The right panel uses a classification based on column densities derived from X-ray observations 
with a dividing column between the groups of 10 21,5 cm -2 . 


The lMerloni et al. f 20140 sample doses not include low redshift very low luminosity AGNs 
and a co uple of notes about this end of the distribution are in order. First, the IR-selected 
sample of Brightman fc Nandra ( 2011bjj shows that the ratio N(type-II)/N(type-I) peaks at 
around L(2-10 keV)~ 10 42 erg s _1 and drops to lower values below this luminosity. Second, 
the discrepancies illustrated in Figure [8] are present in various other studies, including those 
that are based on the measurements of individual covering factors (that are, occasionally, in 
poor agreement between themselves). The sources of confusion have been discussed earlier 
in detail. A major one is, again, the confusion between host galaxy and torus obscuration 
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which is related to the large group of type-Ii sources. 

What is the redshift evolution of obscured and unobscured AGN? Observationally, this 
is very challenging because o f the very small number of observed sources at high redshift. 
Earlier works like Hasineer ( 20081 1 clai m an incre ase in the fraction of obscured sources 
with redshift. This was questioned by Gilli et al. ( 20ld f on grounds of a lack of an ade¬ 


quate K-correction. The I Mer loni ct al.l (12014! ) work, which is perhaps the most complete (a 


statement that will likely be challenged by some X-ray astronomers), is inconclusive about 
this issue regarding low and intermediate luminosity sources. The results for the highest 
luminosity AGNs are more significant with indications for increase in the fraction of ob¬ 
scured sources with redshift amounting to a factor of about 2 between redshifts of 1 and 
3. This Ending is supported by several studies of A GNs in the Chandra Deep Field South 
(CDFS) i Brightman fc IJeda 20121 . Vito et al. 2013^ . 


7.2 Receding and Approaching Tori and Disappearing BLRs 

There were several attempts to explain the dependence of covering factor on source lu¬ 
minosity based on the torus geometry. The idea of a “receding torus”, first suggested by 
Lawrence! 19911 4. involves a very simple torus whose height, h, is independent of the inner 
radius Ri n . Since Ri „ ocLagn 1,/2 , the covering factor is predicted to be proportional to 
Lagn -1 ^ 2 . This ad hoc model makes no attempt to justify the constant height torus. It 
also does not explain the decreasing covering factors deduced from counting X-ray AGNs 
with different obscuring columns. 

Extensive studies based on large AGN samples, some of which have already been men¬ 
tioned, have been used to test the receding torus idea. None of these reproduced, accurately, 
the predicted dependence of fa on Lagn. There were several suggestions to modify the 
original idea. For example, Simpson ! ( 2005 1 measured the fraction of type-I sources by using 
[O hi] A5007 luminosities in SDSS AGNs and constructing LFs for the two types of AGNs. 
The differences between the two LFs were explained to be due to changes in the torus cov¬ 
ering factor with luminosity and were claimed to be removed if h oc Lagn- Unfortunately, 
the work does not take into account the fact that the c onversion between L([0m] A5007) 
and Lagn is luminos it y dep endent ( Stern fc Laorl 2012i l. a fact which was not known at 


lep 

the time. iRic.ci et al.l (1201 31 used the observed EW(Ka) to suggest a somewhat different 


dependence. These attempts, and others, are highly simplified and should not be viewed as 
more than fits of simple analytic functions to the observations. The increased inner radius 
with luminosity is, indeed, a solid and confirmed prediction of all torus models. However, 
other parameters, like fc and h, are not well defined and there is no simple theoretical 
reason to prefer a certain geometry over others. 

The idea of a luminosity dependent sublimation radius has also been used to argue that 
below a certain Lagn, the torus is likely to change its properties, and even disappear. 
This idea deserves the name “ ap p roach ing t ori”. Several of its a s pects were d escribed in 


Laor ( 2003ll . Elitzur fc Shlosman ( 2006l f and Elitzur fc Hoi ( 2009ll . Laor ( 2003h suggested 
to combine the sublimation radius as the outer boundary of the BLR, with the observed 
relationship between BH mass and Lagn- This gives an estimate of the minimal broad line 
velocity (taken here as the full-width at half maximum, FWHM, of the broad H/3 line), 
which is the Keplerian velocity at Ri n - Using current RM results one finds, 


FWHM(H beta ) ~ 1700M g 1/2 L- 1/4 kms- 


(7) 


where, again, I/46=Lagn/10 46 erg s 1 and Mg = Msjt/ 10 8 Mg. For Mg = 1 and 
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Z/46 < 10~ 4 , FWHM(H/3)> 17,000 km s _1 , a line width that is practically undetectable 
in low luminosity AGNs whose continuum near the H/3 line is dominated by stellar light. 
Additional c onsid e ration s show that the line width may be related to LAGN/7/Edd rather 
than Lagn- Laor ( 20031 ') further speculated that there is a natural limit to the maximum 
allowed velocity in the BLR and clouds in this region may not survive close to the BH, 
not just become undetectable. Low luminosity objects that appear to have only narrow 
high ionization lines will be classified as real type-II AGNs and low ionization line sources 
will be classified as type-II LINERs. This results in the increasing fraction of (apparent) 
type-IIs at low luminosities. The argument does not affect the torus covering factor. It 
applies mostly to extremely low lum inosity AGNs in the local universe but not to the real 
type-IIs in the Merloni et al. ( 20141) sample w here Lag n> 10 43 erg s _1 . 

Disappearing BLRs have been discussed by Nicastro! ( 2000 1 who suggested that gas out¬ 
flow from the surface of the disk, outside the regio n dominated by radi ation pressure, can 
lead to disap pearing BLRs at very low L AGN/LEdd- Elitzur fc Ho 1 2009h followed the ideas 
developed in Elitzur fc Shlosmanl ( 2006l l and considered the clumpy dust-free BLR, and 
the clumpy dusty torus, as one entity made of clouds lifted along magnetic field lines from 
the surface of the disk. A limiting radial column density of ~ 10 22 cm~ 2 , was assumed 
to be a requirement for producing observable broad emission lines from the clumps. The 
work shows that at low mass accretion rate, corresponding to Lagn^ 10 40 , BLR clouds, 
and large dusty clumps, no longer exist. This luminosity, which varies with the black-hole 
mass, and the corresponding I/AGN/LEdd, are i n the right ball park for many LINERs but is 
several orders of magnitude smaller than the typical luminosity of many real type-II AGNs 
( Brightman fc Nandral 201 lbl . Marinucci et all 20121 . Merloni et al. 20hdl . Despite the dis¬ 
crepancy, the general idea can perhaps be extended to include more luminous sources with 
no detectable broad emission lines. 


7.3 Mergers Evolution and Unification 


Detailed morphological studies of host galaxies of AGN s, up to redshift ~ 2, reveal that 
about ~ 80% of them do not show signs of interaction (e.g. Cisternas et al . 201ll. ISchawinski et al 


20121 1. Such studies depend on signal-to-noise, spatial resolution, and most importantly red- 
shift. They should be interpreted as strong evidence against major mergers and perhaps 
not so strong evidence regarding the final stages of a minor merger where the mass ratio of 
the two hosts is large. The finding, which contradicts earlier claims about the association of 
most AGN activity with galaxy mergers, is evidence for the importance of secular processes 
in transporting cold gas from the disk to the vicinity of the BH. The studies show also that 
a large fraction of the most luminous AGNs at high redshift are associated with mergers. 
Indeed, merger induced nuclear activity is likely to be more common at very high redshift, 
when the density of the universe was higher and much of cold gas was not yet converted 
to stars. Can simple unification based on relatively simple shaped tori, and well ordered 
motion of cold gas, be applied to the complex geometry of such mergers? How early in the 
history of an AGN triggered by a major merger can such an ordered central structure form? 
Can merger driven AGNs be classifi ed like the secularly evolving systems if their BLR is 
obscured by galactic scale dust (e.g. Imanishi. Maiolino fc Nakagawa 201 0)1. and what are 
the implications to the NIR-MIR SEDs of such systems? 

Major mergers of gas rich systems in the early universe have been suggested to result 
in powerful dusty SF systems (dust-obscu red galaxies, DO Gs) su rrounding a buried AGN 


whose power is increasing with time (see Narayanan et alJhoiol . and references therein). 
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As the AGN power increases, the surrounding gas is blown away and an unobscured AGN 
is revealed. The most luminous AGNs during the era of peak AGN activity, at z=2- 
3, have been suggested to be objects of this type. The geometry of the obscuring dust 
in this case must be very different from the simple tori assumed in non-merging galax¬ 
ies. Unfortunately, numerical merger simulations do not yet have sufficient resolution to 
model the geometry of such structures and high resolution imaging is limited to nearby 
systems (see however Schawinski et al. ( 2012h ~l. Thus, SED observations and model fit¬ 
ting are the best tools to investigate these scenarios, in particular in objects where both 
optical and IR information is available (e.g. SDSS and WISE see Weedman et al. 20121 ; 


Vardanyan. Weedman fe Sargsvanl 2014). 


An empirical way to study such scenarios is to examine objects with very large Ldust /Lagn 
and check whether this is caused by UV luminosity which is unusually low (UV extinction 
in a DOG) or MIR luminosity w hich is unusually high (large coveri n g fact or and a clear 
LOS to the center). The work of Vardanyan . Weedman fc Sargsvan ( 20141 '! suggests that 
the latter is correct and shows that the objects with the smallest Ldust /Lagn are those with 
the largest L(UV). Such studies show no evidence of the DOG-related scenario in z=2-3 
AGNs (see however Glikman et al. ( 2012h for evidence to the contrary). 

There are more detailed studies of merging systems at low and intermediate redshifts, 
especially ultra luminous infrared galaxies (ULIRGs) hosting active BHs. A low red- 
shift sample of this type i s the QUEST (Qua sar/ULIRG Evolutionary Study) sample (see 


Veilleux et al. 20061 . 20091 . Veilleux et al. 20091 . and references therein). Some of the sources 


in the sample contain low luminosity AGNs, including LINERs, and are clearly SF domi¬ 
nated systems. Others contain high luminosity AGNs, mostly PG Quasars.. HST imaging 
of QUEST sources indicate that the type of AGN depends on the evolutionary stage of 
the host, and the BH activity is most likely triggered by the merger. In particular, there 
is evidence that a merger progresses from type-II at the earlier stages, to type-I closer 
to its end. Type-I AGNs become more and more com mon among fully-coalesced merger 
remnants with elliptical-like hosts I Veilleux et al. 2006). They also become more common 


with increasing Ldust- This may or may not be related to the general trend of decreas¬ 
ing obscuration with increasing L(2-10 keV) discussed for non-merging systems. Many of 
the QUEST sources have been observed, spectroscopically, b y Spitzer. Interestingly, their 
comparison with sou rces that show no indication of mergers ( Mor. Netzer fc Elitzur 20091 : 
Mor fe Netzer 20121 '! shows very similar SEDs and a similar trend of decreasing covering 


factor with Lagn- 

Can major mergers at high redshift produce galactic-scale obscuration that results in torus 
like SED and absorption? Can heating by the AGN result in cold (30-50 K) dust emission 
that mimics SF regio ns? There are no clear observational answers to these questions but 
recent calculations bv ISchneider et all (2014) are already giving some clues. 

A complex geometry that can give rise to both MIR and FIR emission, is a “torus within 
a torus” or “an obscurer within an obscurer”. In this scenario, a large covering factor 
obscurer, just outside the dust sublimation radius, results in intense NIR-MIR emission. 
This obscurer is clumpy and has a simple axisymmetric geometry. Part of the AGN radiation 
escapes through its ionization cones and part leaks out through the clumps. A much larger 
obscurer, with an undefined geometry that does not block the central dust-free cone, is 
present at a distance of a few kpc. The covering factor of the inner obscurer is fc( in) and 
the one of the outer obscurer fc{ out). It is easy to think of several scenarios that result in 
different relative L(FIR)/L(MIR)/Lagn, depending on how much radiation of the central 
source leaks through the small obscurer near the center. It is also possible that some of the 
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MIR radiation of the central obscurer is absorbed by the further dust, if its optical depth is 


large enough. A certain variant of this complex geometry was modeled bv ISchneider et al. 
(2014) who found that such a geometry with /c(in)~ 0.7, gives a reasonable fit to the 
observed SED of a specific z=6.4 AGN. In that particular model, about half of the FIR. 
luminosity is AGN-related and half SF-related. 

An obvious observational test would be to obtain polarimetry maps of such systems since 
this is a most effective way to check whether the distribution of the position angles is 
centrosymmetric (i.e. scattering from a central point) or not (see Kishimoto 19991 ). 

In summary, it is hard to imagine that major mergers, especially in the early universe, 
can lead to the highly ordered geometry, and the torus-like structure, suggested to explain 
the observations of lower luminosity, secularly evolving systems. Detailed SED studies of 
high redshift mergers can help to resolve this issue. 


8 SUMMARY: PRESENT AND FUTURE UNIFICATION 

The unification scheme presented in this review is different in some ways, and similar 
in others, to the scheme presented in the late 1980 that had three fundamental pillars: 
orientation, covering factor and luminosity. The field is driven by courageous attempts to 
construct simple models for complex situations, and attempts to understand sever selection 
effects typical of large samples. The most important developments of the last decade or so 
can be summarized as follows: 

1. Most AGN tori are probably axisymmetric and clumpy. Unfortunately, todays phe¬ 
nomenological clumpy torus models leave much to be desired and their parameters are 
hard to constrain. Hydrodynamical calculations that follow the inflow of gas from the 
galaxy, are more promising. They naturally produce a clumpy thick central structure 
and demonstrate the importance of SN explosions, and radiation pressure feedback 
by the central AGN. The simulations do not yet include all relevant processes and are 
generally limited to AGNs with small BHs. 

2. Disk-wind models of different types seem to be able to explain several fundamental 
properties of the central obscurer, like the formation of dust-free and dusty clumps 
at roughly the desired locations. So far, there is no clear preference for magnetic or 
radiation pressure driven winds and there is a general problem, in all models, to form 
and maintain large column density clumps. 

3. The central issues of torus covering factor and anisotropy are not yet fully understood. 
The general trend of decreasing fc with increasing Lagn is confirmed but its details 
not well understood. There are substantial disagreements between the various meth¬ 
ods used to study this relationship, like the ones based on Ldust/LAGN, on EW(Ka), 
or the division of the population into groups based on X-ray absorbing columns. The 
suggestion of a receding torus is the only idea driving this field. However, the idea is 
neither physically sound nor well understood. 

4. Study of ionization cones and IR interferometry are two major developments in this 
field. Bicone observations provide the most reliable way to define the torus bound¬ 
aries and link the emission line properties to the kinematic properties of the NLR 
gas. Long baseline interferometry is the way of the future. It has already provided 
global intriguing results, like an unexpected run of dust temperature with radius and 
model-images that indicate deviations from a simple torus-shaped obscurer. They 
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also stress the need to reconsider the role of MIR emitting polar outflows. Future 
phase information will likely provide real images of central AGN tori. 

5. Part of the confusion in the present unification scheme is the presence of several 
sub-groups that may not belong in it in the first place. In particular, “real type- 
II AGNs”, that are not really understood, confuse several issues like the fraction of 
obscured and unobscured low luminosity AGNs. Another group that shows a mixed 
bag of properties are type-Ii AGNs. Some of these sources are heavily extinguished 
by dust in the host galaxy, others are type-Is “caught” during a time of minimum 
AGN activity. 

6. Unification does not seem to stand the test of BH evolution. In particular, merging 
systems likely behave in a different way from AGNs in secularly evolving hosts. Based 
on todays observations, there is no reason to believe that accreting BHs in major 
mergers, or merger remnants, partake in the simple, old unification scheme. 

Returning to the central theme of unification, it is clear that the old scheme requires several 
major revisions. The more important ones are the recognition that for a given luminosity 
there is a large range of covering factors, the realization that X-ray absorption can be caused 
by two and possibly more components, and the identification of several sub-groups of AGNs 
that do not belong in this scheme in the first place. A revised scheme would require better 
ways to isolate and remove cases where obscuration has nothing to do with the nucleus, 
and improved X-ray techniques to clarify the importance of X-ray absorption by dust-free 
gas. Perhaps the most crucial aspect, which is observationally challenging, is the separation 
of torus absorption from the presence of broad emission lines. Putting the absorption of 
the non-stellar continuum in the center of the scheme, and minimizing the use of broad 
emission lines as a major characteristics, are good ways to proceed. 
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